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FOREWORD 

This  report  was  prepared  by  1st  Lt.  George  R.  Frost  and 
Dr.  Arthur  J.  Wennerstrom  of  the  Fluid  Dynamics  Facilities 
Research  Laboratory,  Aerospace  Research  Laboratories,  Wright- 
Fatterson  Air  Force  Base,  Ohio,  and  by  Hr.  Richard  M.  Hearsey 
of  the  University  of  Dayton  Research  Institute. 

The  report  presents  results  from  a  portion  of  the  effort 
of  the  Fluid  Machinery  Research  Group,  supervised  by  Dr.  Arthur 
J.  Wennerstrom  and  was  conducted  under  Work  Unit  09  of  Project 
7065,  "Aerospace  Simulation  Techniques  Research"  under  the 
overall  direction  of  Mr.  Elmer  G.  Johnson. 


ABSTRACT 


This  report  is  a  revision  and  extension  of  ARL  70-00A6. 

It  describes  the  analysis  in,  and  the  use  of,  a  computer 
program  which  has  been  developed  for  use  in  the  design  of  axial 
compressor  airfoils  suitable  for  operation  at  high  subsonic  and 
supersonic  Mach  numbers.  Four  rather  versatile  camber  line 
shapes  and  two  thickness  distributions  are  mathematically 
derived.  These  camber  lines  provide  the  capability  of  defining 
a  wide  variety  of  blades,  from  those  of  continuously  positive 
camber  to  the  so-called  "S-blades",  including  many  of  the 
intermediate  possibilities.  A  method  is  presented  whereby  the 
airfoils  are  specified  on  arbitrary  axisymmetric  streamsurfaces 
and  then  accurately  redetermined  in  Cartesian  coordinates  on 
planes  normal  to  the  stacking  axis.  The  program  determines  the 
coordinates,  calculates  section  properties,  and  provides 
precision  plots  of  the  specified  blade  sections  and  the  result¬ 
ing  blade  projections  on  the  above-mentioned  planes. 
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SECTION  I 


INTRODUCTION 


The  Fluid  Dynamics  Facilities  Research  Laboratory  of  the 
Aerospace  Research  Laboratories  (ARL)  is  engaged  in  a  program 
to  find  means  of  designing  axial  compressors  capable  of  efficient 
performance  at  Mach  numbers  and  pressure  ratios  beyond  the 
current  state-of-the-art.  In  conjunction  with  this  effort,  it 
was  determined  that  general  improvements  could  be  made  over 
some  of  the  blade-design  techniques  now  in  popular  use  for  high 
relative  Mach  numbers. 

It  is  desirable  when  designing  a  compressor  blade  to  create 
blade  sections  upon  arbitrary  surfaces  of  revolution  that 
correspond  to  streamsurfaces  indicated  by  the  compressor  design 
calculations.  However,  for  manufacturing  purposes,  it  is 
desirable  to  have  blade  surface  coordinates  for  sections 
through  the  blade  that  are  plane  and  perpendicular  to  the 
stacking  axis.  A  method  is  presented  whereby  the  blade  is 
designed  by  creating  sections  on  the  streamsurfaces,  and  then 
determining  coordinates  for  plane  sections  through  the  resultant 
blade. 

Several  objectives  guide  the  aerodynamic  design  of  any 
airfoil,  whether  for  a  compressor  or  any  other  purpose.  First, 
and  usually  foremost,  is  the  desire  to  keep  the  rate  of 
diffusion  as  low  as  possible  on  the  airfoil  surfaces.  While 
this  objective  does  not  change  with  Mach  number,  the  shape  of 
airfoil  required  to  accomplish  this  can  change  significantly 
with  Mach  number.  At  supersonic  Mach  numbers  there  is  generally 
little  camber,  and  sometimes  negative  camber,  in  the  leading 
edge  region  of  a  compressor  blade.  At  these  Mach  numbers, 
minimizing  diffusion  rates  means  minimizing  the  suction  surface 
supersonic  expansion  upstream  of  the  inevitable  passage  shock 
or  even  producing  a  limited  amount  of  supersonic  diffusion 
upstream  of  the  shock  location.  There  are  two  principal  methods 
of  obtaining  supersonic  diffusion.  When  using  blades  having 
zero  or  positive  camber  near  the  leading  edge,  the  diffusion 
may  be  obtained  by  a  convergence  of  the  annulus  walls.  If  this 
is  impractical,  for  example,  in  a  high-aspect-ratio  supersonic 
fan  rotor,  the  same  result  can  be  accomplished  by  negative 
camber  in  the  leading  edge  region  of  the  airfoil.  Blades  which 
have  negative  camber  followed  by  a  region  of  positive  camber 
have  been  called  "S-blades". 

A  second  objective  is  to  define  the  airfoil  in  such  a 
manner  that  its  shape  can  be  continuously  varied  to  accommodate 
varying  conditions  along  the  blade  span,  thus  maintaining 
resonable  diffusion  rates  along  each  blade  element.  For 
example,  it  would  be  possible  to  define  a  type  of  profile 
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having  a  very  attractive  diffusion  distribution  for  a  hub  section, 
but  which  could  not  be  adapted  to  tip  conditions,  or  vice  versa. 
This  would  obviously  be  undesirable.  It  would  generally  be 
impractical  to  attempt  to  blend  airfoil  sections  of  more  than 
one  type  along  the  span  of  a  blade.  A  third  objective  for  a 
desirable  compressor  blade  is  that  the  parameters  determining 
its  aerodynamic  characteristics  should  be  variable  over  a  wide 
range  within  which  there  is  little  danger  of  producing  shapes 
that  are  mechanically  undesirable. 

This  report  is  a  revision  and  extension  of  an  earlier 
report  (Reference  1)  which  presented  two  blade  profiles  which 
satisfy  all  the  above  requirements.  They  are  composed  of  a 
common  thickness  distribution  applied  to  two  different  camber 
lines:  one  a  fourth-order  polynomial,  the  other  defined  by  two 
exponential  functions.  The  thickness  distribution  consists  of 
two  third-order  polynomials:  one  between  the  leading  edge  and 
the  point  of  maximum  thickness,  and  one  from  there  to  the 
trailing  edge. 

In  addition  to  these  camber  lines,  two  other  camber  lines 
have  been  incorporated  into  the  associated  computer  program. 

One  of  these  is  the  familiar  circular-arc  camber  line,  upon 
which  is  placed  the  double-circular-arc  thickness  distribution. 

The  second  camber  line  may  be  composed  of  any  pair  of  circular 
arcs,  either  of  which  may  degenerate  to  a  straight  line, 
allowing  for  a  wide  range  of  design  flexibility.  This  latter 
camber  line,  referred  to  as  the  multiple-circular-arc  camber 
line  in  this  report,  uses  the  common  thickness  distribution  of 
the  original  blade  profiles  and,  as  a  special  case,  allows 
this  thickness  distribution  to  be  applied  to  a  simple  circular- 
arc  camber  line . 

The  various  camber  lines  and  thickness  distributions  are 
described  in  Sections  II  and  III. 

In  order  to  specify  a  blade  for  manufacture  it  is  desirable 
to  have  coordinates  for  plane  sections  through  the  blade.  For 
maximum  accuracy,  blade  profiles  should  be  designed  on  the 
arbitrary  surfaces  of  revolution  that  typical  current  design 
procedures  indicate.  The  method  of  satisfying  these  two 
requirements  consists  first  of  designing  blade  profiles  on 
surfaces  of  revolution.  Then  a  thin,  uni formly-f lexible  beam 
is  passed  (mathematically)  through  like  points  on  each  profile 
and  the  intersection  of  the  beam  with  the  "manufacturing 
planes"  yields  the  coordinates  of  the  "manufacturing  sections," 

The  original  report  and  computer  program  presented  Equation  (Ji?) 
Reference  1,  incorrectly  which  resulted  in  an  improper 
determination  of  the  Cartesian  coordinates  of  the  specified 
streamsurface  blade  sections  for  any  non-cy lindri cal  stream- 
surface.  This  led  to  associated  errors  in  the  determination  of 
the  manufacturing  sections.  (The  magnitude  of  this  error 


increases  with  streamsurface  slope  angle  and  generally  would 
not  result  in  a  blade  leading  or  trailing  edge  angle  error  of 
more  than  about  one  degree.)  The  version  cf  the  program  present¬ 
ed  in  this  report  incorporates  the  correct  analysis  as  described 
in  Section  IV. 

The  fifth  section  of  this  report  describes  calculations  to 
determine  various  properties  of  both  the  streamsurface  blade 
profiles  and  the  manufacturing  sections.  Several  blade 
characteristics  required  for  aerodynamic  analysis  of  a  resulting 
blade  row  are  presented  in  Section  VI. 

In  addition  to  the  new  features  and  the  analysis  correction 
already  mentioned,  the  computer  program  (presented  in  Sections 
VII  through  17)  has  been  modified  in  several  other  respects 
since  the  original  report  was  prepared.  The  program  now 
accommodates  curvilinear  computing  stations,  as  compared  to 
the  strictly  linear  stations  previously  required.  An  additional 
section  property,  the  torsional  constant,  is  calculated  for 
the  manufacturing  sections  to  facilitate  structural  analysis  of 
the  blade.  Finally,  a  blade  may  now  be  stacked  also  at  either 
the  leading  or  trailing  edge,  in  addition  to  the  previous 
capability  of  stacking  at,  or  offset  from,  the  section  centroid. 

The  seventh  section  of  this  report  presents  a  general 
user's  guide  for  the  computer  program,  and  Section  VIII  presents 
several  examples  for  the  use  of  the  program.  The  actual  coding 
and  a  synopsis  of  the  essential  steps  in  the  program  logic, 
as  well  as  information  which  may  facilitate  its  implementation 
or,  a  computing  system,  are  presented  in  Section  IX. 
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on  a  computing  system,  are  presented  in  Section  IX. 


SECTION  II 


THE  SECTION  CAMBER  LINE 


To  obtain  a  camber  line  varying  smoothly  in  curvature  and 
satisfying  the  requirements  of  a  compressor  blade,  it  would  be 
convenient  to  describe  it  mathematically  in  equation  form.  The 
type  and  order  of  such  an  equation  will  in  turn  be  influenced 
by  the  boundary  conditions  which  are  imposed.  The  equation 
must  be  of  high  enough  order  to  accommodate  all  of  the 
necessary  boundary  conditions.  Its  order  should  be  no  higher 
than  absolutely  necessary,  however,  to  minimize  the  number  of 
singularities  which  can  occur  within  its  useful  working  range. 

In  the  following  discussion,  we  shall  consider  an  equation  of 
the  form  y  =  f(x)  where  the  variable  x  is  nondimensional  and 
varies  from  0  at  the  blade  leading  edge  to  1.0  at  the  blade 
trailing  edge.  Furthermore,  instead  of  defining  the  x-axis  as 
the  chord  line,  which  is  usual,  x  has  been  defined  as  the  axial 
direction  in  the  blade-element  plane.  This  definition  of  the 
x-axis  was  chosen  because  it  was  found  to  lead  to  the  simplest 
treatment  of  boundary  conditions. 

Three  boundary  conditions,  or  their  equivalent,  are 
absolutely  necessary.  One  condition  fixes  a  point  on  the  line 
with  respect  to  the  coordinate  system,  another  defines  the 
slope  at  the  leading  edge,  and  the  third  defines  the  slope  at 
the  trailing  edge.  These  three  alone  are  sufficient  to  define 
a  unique  circular-arc  camber  line  between  the  edge  points. 
However,  these  boundary  conditions  are  insufficient  to  define 
a  camber  line  for  a  blade  which  will  meet  the  special  require¬ 
ments  of  efficient  operation  at  relatively  high  Mach  numbers. 

A  characteristic  usually  desired  for  supersonic  compressor 
blades  is  very  little  camber  in  the  leading  portion  of  the 
blade.  Depending  upon  the  overall  aerodynamic  design,  the 
desired  leading-edge  camber  may  be  negative,  zero,  or  slightly 
positive.  A  convenient  way  of  controlling  this  is  to  impose 
a  boundary  condition  on  the  second  derivative  of  the  camber 
line  at  the  leading  edge.  To  put  it  in  dimensionless  form,  we 
specify  the  ratio  between  second  derivatives  at  the  leading 
edge  and  the  point  where  the  absolute  value  of  the  second 
derivative  is  a  maximum.  A  fifth  condition  was  felt  advisable 
for  an  efficient  compressor  blade.  Namely,  the  first  four 
conditions  used  alone  tend  to  produce  a  blade  whose  curvature 
is  highest  at  the  trailing  edge.  This  could  result  in  large 
deviation  angles  and  correspondingly  high  losses.  Therefore, 
a  condition  was  imposed  on  the  second  derivative  at  the  trailing 
edge  whereby  its  value  could  be  specified  as  something  less 
than  the  maximum.  This  boundary  condition  was  phrased  in 
exactly  the  same  manner  as  at  the  leading  edge.  The  resulting 
five  conditions  may  be  written 


.’t 


At  x  =  0:  y  =  0 

y‘  =  tan  ctj 

y"  =  P  (y" ) 

At  x  =  1 :  y '  =  tan  a2 

y"  =  Q  (y") 


max 


(1) 


(2) 


max 


where  prime  (')  and  double  prime  (")  refer  to  first  and  second 
derivatives  respectively. 

1.  POLYNOMIAL  CAMBER  LINE 

The  five  stated  boundary  conditions  could  be  applied  to  a 
great  number  of  equations  to  produce  a  camber  line.  A  simple 
fourth  order  polynomial  was  chosen  for  the  blade  having 
primarily  positive  camber  because  it  was  convenient,  and  proved 
quite  satisfactory.  Since  the  second  derivative  of  a  fourth 
order  polynomial  is  obviously  second  order  and  two  of  the 
boundary  conditions  apply  to  the  second  derivative,  it  was 
convenient  to  start  by  phrasing  this  in  the  form  of  a  parabola. 
In  standard  form  this  can  be  written 

(x  -  h)2  =  4a  (y"  -  k)  (3) 

or  l 

y"=— (x-h)2+k 

where  h  is  the  position  on  the  x-axis  where  the  second 
derivative  is  a  maximum  (or  minimum)  and  k  is  the  value  which 
the  second  derivative  has  at  this  point.  Integrating  twice, 
one  obtains 

y'=Y^-(x-h)3  +  kx  +  b 

y  =  -1-  (x  -  h)4  +  -  xz  +  bx  +  c  (5) 

y  48a  v  2 

Application  of  the  five  boundary  conditions  to  these  equations 
produces  the  following  result 

At  x  =  0,  y  =  0  (g) 


c 


h4 

48a 
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If  1-p 


The  values  of  the  input  parameters  at  which  the  equations 
degenerate  are  easily  ascertained  by  examining  the  derived 
constants.  The  values  of  the  two  blade  angles  may  be  anything 
between,  and  less  than,  ±90  degrees,  measured  from  the  axial 
direction.  Mote  that  when  the  two  blade  angles  are  the  same, 
this  fourth-order  camber  line  becomes  a  straight  line, 
regardless  of  what  values  P  and  Q  may  have.  The  parameter  Q 
may  be  given  any  value  between  zero  and  1.0.  In  practice,  a 
value  of  0.5  has  proven  very  satisfactory  for  most  applications. 
Singularities  associated  with  the  parameter  P  are  partially 
dependent  upon  Q.  At  P  =  1.0,  a  singularity  exists  independently 
of  Q.  However,  if  Q  =  1.0,  the  equations  also  degenerate  at 
P  =  -2.0  and  if  Q  =  0,  the  equations  degenerate  at  P  =  -3-0. 

In  general,  P  should  be  greater  than  -2.0  and  less  than  1.0. 

For  most  applications,  a  value  of  P  =  0  is  quite  satisfactory. 
Toward  the  hub  of  a  compressor,  it  may  be  desirable  to  give  P 
small  positive  values.  Toward  the  tip,  one  might  wish  to  use 
slightly  negative  values  of  P.  However,  this  particular 
camber  line  does  not  produce  an  attractive  "S-blade,"  especially 
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if  the  overall  section  camber  is  nearly  aero.  Negative  values 
of  P  are  only  useful  for  delaying  the  onset  of  positive  curvature 
by  producing  very  slight  negative  cambers  at  the  leading  edge. 

Figures  1  through  5  illustrate  the  general  characteristics 
of  this  camber  line.  In  Figure  1,  the  distribution  of  the  second 
derivative  of  the  camber  line  is  shown  for  various  values  of  P 
and  blade  angles  characteristic  of  mid  span.  Figure  2  shows 
the  actual  camber  lines  which  correspond  to  the  curves  shown  in 
Figure  1.  Figure  3  is  similar  to  Figure  2  but  shows  the 
influence  of  P  on  a  camber  line  more  typical  of  a  compressor 
hub  section.  Figure  ^  illustrates  the  influence  of  Q  for  a 
camber  line  typical  of  mid  span  with  P  =  zero.  Finally, 

Figure  5  shows  the  different  camber  lines  which  result  from 
varying  the  blade  outlet  angle  with  all  other  parameters  fixed 
at  typical  values. 

2.  EXPONENTIAL  CAMBER  LINE 

The  exponential  camber  line  is  the  result  of  efforts  to 
find  a  mathematical  expression  capable  of  defining  a  satisfactory 
"S-blade."  An  "S-blade"  should  satisfy  all  of  the  basic 
criteria  which  led  to  the  preceding  polynomial  camber  line  and, 
in  addition,  several  conditions  unique  to  the  !,S"  configuration. 
Namely,  it  should  be  possible  to  locate  the  inflection  point 
anywhere  on  the  camber  line,  the  blade  angle  at  that  point 
should  be  independently  specifiable,  and  the  method  must  work 
for  overall  blade  cambers  which  are  positive,  zero,  or  negative. 
Furthermore,  the  transition  across  the  inflection  point  should 
occur  smoothly  but  rapidly  so  that  a  long  straight  region  is 
avoided  in  the  middle  of  the  blade.  It  is  also  highly  desirable 
that  the  equations  allow  a  smooth  transition  from  the  "S"  config¬ 
uration  to  a  more  conventional  one  such  as  characterized  by 
the  polynomial  camber  line.  To  the  list  of  five  boundary 
conditions  presented  earlier,  one  is  added: 

At  x  =  s ,  y 1  =  tan  as  (11) 


Numerous  equations,  including  polynomials,  were  examined 
with  respect  to  the  six  boundary  conditions.  It  soon  became 
evident  that  more  satisfactory  results  could  be  obtained  by 
defining  the  camber  line  with  two  equations,  one  either  side 
of  the  inflection  point.  As  with  the  previous  blade,  it  was 
easiest  to  choose  an  equation  by  examining  equations  for  the 
second  derivative  of  the  camber  line  and  then  integrating  twice 
to  obtain  the  coordinate  equation.  The  equation  selected  for 
both  portions  of  the  camber  line  is 


y 


II 


b  ( x 


s) 


a( x-s ) 

e 


(12) 
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One  set  of  constants  applies  from  0  to  s  and  another  set  applies 
from  s  to  1.0.  The  maximum  value  of  y"  occurs  where  y’*'  =  0, 
which  is  at 

x  =  s  -  -  (13) 

a 


Note  that  the  constant  a  is  positive  ahead  of  the  inflection 
point  s,  and  is  negative  after  it.  The  maximum  (or  minimum) 
value  of  the  camber  line  second  derivative  for  each  camber  line 
segment  is 


y 


II 

max 


_b_ 

ae 


(14) 


Integrating  the  second-derivative  equation  twice,  one  obtains 


y ' 


4  ea<x-s)  [a(x-s)  -  1]  +  c 

a  2 


(15) 


b_  ga(x-s)  [a(x_sj  -  2]  +  c  (x-s)  +  d 


(16) 


This  equation  contains  four  arbitrary  constants,  and  there  are 
two  sets  to  be  determined  as  described  above.  Because  there  are 
different  equations  describing  the  two  portions  of  the  camber 
line,  a  seventh  condition  is  added  to  the  six  already  stated. 

At  the  inflection  point,  the  coordinates  of  the  camber  line  are 
the  same  on  both  lines . 

Applying  the  four  conditions  appropriate  for  the  leading 
portion  of  the  camber  line  produces  the  following  result: 

At  x  =  0,  y  =  0 

b , 

d1  =  ( a !  s  +  2) —  e"ais  (17) 

ai 
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At  x  =  0,  y'  =  tan  ctj 

Cj  =  tan  a!  +  UjS  +  1 )  — e"ais 

ai 


At  x  =  0, 

y"  =  P  (y"  ) 

max' 

n  ( 1-a, s ) 

P  =  a, se  1 

(18) 

At 

x  =  s ,  y 1  =  tan  ag 

(19) 

b, 

c,  =  tan  a  +  — r 

1  s  a. 

Equating  the  two  expressions  for  ct  yields 

a,2(tan  a.  -  tan  a  )  (20) 

b  =  - l - 

1  -  ( a j s  +  l)e"ais 

Applying  the  four  conditions  appropriate  for  the  rearward 
portion  of  the  camber  line  produces  the  following  result: 

At  x  =  s,  y '  =  tan  a 
J  s 

c2  =  tan  a  +  — \ 

2  s  a, 2 


At  x  =  s,  y  =  y  obtained  from  Equation  (16)  for  the  leading 
portion  of  the  camber  line. 


At  x  =  1  ,  y '  =  tan  a2 


c 2  =  tan  a2  -  ea2(1"s)  [a2 ( 1 -s )  -  1]  (22) 

a2 
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(23) 


Equating  the  two  expressions  for  c2  yields 

a,2(tan  a,  -  tan  a  ) 
f _ _ _ s _ 

2  1  +  [a2 ( 1 -s)  -  1]  ea2(1_s) 


At  x  -  1 ,  y"  =  Q  (y "  ) 

*  J  y  w  max' 


Q  = 


(s-l)a2e1+a2(1“s) 


(24) 


The  equations  obtained  to  determine  a,  and  a2  (Equations 
(18)  and  (24))are  implicit,  and  require  iterative  solution. 

In  each  case,  there  are  two  solutions  for  the  constant  (assuming 
P  and  Q  to  be  less  than  unity).  Because  it  is  desired  that  the 
magnitude  of  the  second  derivative  reach  its  maximum  value 
between  the  inflection  point  and  leading  or  trailing  edge,  as 
appropriate,  it  follows  that 


and 


> 


2 

s 


< 


1 

s-1 


(25) 


Figures  6  through  10  illustrate  the  general  characteristic 
of  this  camber  line.  Figure  6  shows  the  distribution  of  the 
second  derivative  of  the  camber  line  for  various  values  of  P  and 
Q.  The  other  parameters  that  are  held  constant  describe  a 
blade  perhaps  typical  of  mid-span  conditions  for  a  high-speed 
machine.  Figure  7  shows  the  corresponding  camber  line. 

Figure  8  shows  camber  lines  for  a  section  having  zero  overall 
camber  and  hence  equal  negative  and  positive  cambers.  The 
effect  of  varying  the  parameters  P  and  Q  is  shown.  Figure  9 
shows  camber  lines  for  a  similar  section,  except  that  the  effect 
of  varying  the  location  of  the  inflection  point  S  is  illustrated. 
Camber  lines  for  three  sections  having  no  negative  camber  are 
shown  in  Figure  10.  As  the  "inflection  point"  is  shifted 
rearwards  from  the  leading  edge,  the  camber  line  changes  from 
one  of  continuous  positive  camber  to  that  of  a  "J-blade." 

3.  CIRCULAR-ARC  CAMBER  LINE 

The  circular-arc  camber  line  has  been  considered  because 
of  the  wide  application  of  the  double-circular-arc  airfoil  in 
turbo-machinery  of  the  more  conventional  variety  and  because 
of  its  potential  use,  particularly  with  respect  to  stators,  in 
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high  performance  compressor  stages.  Of  the  original  five  boundary 
conditions  of  Equations  (1)  and  (2),  those  involving  y"  are 
superfluous  to  the  definition  of  a  unique  circular-arc  camber 
line. 


The  equation  of  this  camber  line  is  of  the  form 

(x  -  xq ) 2  +  (y  -  y0)2  =  R2  (26) 


where  (x0,  y0)  is  the  center  of,  and  R  the  radius  of,  the  circle 
of  which  the  desired  arc  is  a  portion. 

From  the  boundary  condition  at  x  =  1  and  the  knowledge  that 
the  slope  of  a  radial  line  at  any  point  must  equal  the  negative 
reciprocal  of  the  slope  of  the  tangent  to  the  arc  at  that  point, 
one  obtains 


tan  a2 


(27) 


where  y2  is  the  y  coordinate  of  the  arc  at  x  =  1.  Solving 
for  x0, 

xQ  =  (y2  -  yQ)  tan  a2  +  1  (28> 


At  the  point 


(1, 


y2),  Equation 

( y 2  -  y0)2  = 


(26)  may  be  arranged  to  yield 

-  V  +  2*o  - 1 


(29) 


From  the  boundary  condition  which  fixes  the  leading  edge 
of  the  arc  at  the  origin  of  the  coordinate  system,  and  Equation 
(26),  one  obtains 


y 


o 


2  _ 


(30) 


Substituting  Equations 


iy2  -  y0)2 


(28)  and  (30)  into  Equation  (29), 
=  yQz  +  2  (y 2  -  yQ)  tan  a2  +  1 


(3D 
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Expanding  the  left  hand  side  of  Equation  (31)  and  solving 
for  y0 , 


yn  = 


y22  -  2 y2  tan  a2  -  1 
2(y2  ~  tan  <x2) 


Using  this  expression  for  yQ  in  Equation  (28), 


xo  = 


(y22  -  1)  tan  <*2  +  2y. 
2(yz  -  tan  a2) 


and  substituting  Equations  (32)  and  (33)  into  (30)  and  solving 
for  R2 , 


R2  = 


( y 2 2  +  I)2  sec2  ct2 
4(y2  -  tan  a2)2 


The  above  equations  (Equations  (32),  (33)  and  (3*0)  give 
the  necessary  camber  line  constants  in  terms  of  the  as  yet 
undetermined  y2 .  The  third  boundary  condition,  fixing  the 
slope  at  the  origin,  is  employed  to  express  y0  in  terms  of 
known  quantities.  From  this  boundary  condition, 

xo  =  "yo  tan  ai  (: 


Combining  Equations  (28)  and  (33)  and  solving  for  y2 , 

yQ  (tan  a2  -  tan  ax)  -  1 

y  z  7 

tan  a„ 


(36) 


Substituting  this  expression  for  y2  into  Equation  (32) 
yields,  after  considerable  simplification, 


yo  = 


-  tan  aj  sec2  a2  ±  tan  a2  sec  <Xj  sec  a2 


tan2  Oj  -  tan2  a2 
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The  proper  sign  of  the  second  term  (y)  of  Equation  (37) 
is  determined  by  the  realization  that  J  y  j  must  approach 
infinity  as  a2  approaches  ctx .  Therefore, 

-  (tan  ctj  sec2  a2  +  tan  ax  sec  sec  a2)  (38) 

y  =  - - 

tan2  cij  -  tan2  a2 


Substitution  of  Equation  (38)  into  Equation  (35)  determines 
xQ  in  terms  of  known  quantities.  Substitution  of  Equation  (38) 
into  (36)  determines  y2  in  terms  of  the  known  quantities  0^ 
and  a2 ,  which,  when  used  in  Equation  (3*0,  yields  an  1 
expression  for  P?  in  terms  of  a,  and  a2  .  This  completes  the 
determination  of  the  three  constants  (x  ,  y  ,  and  R2)  of  the 
circular  arc  camber  line. 

MULTIPLE-CIRCULAR-ARC  CAMBER  LINE 

The  multiple-circular-arc  camber  line  incorporates  and 
extends  the  circular-arc  definition  scheme  as  previously 
described  to  include  the  potential  for  defining  a  camber  line 
composed  of  two  segments,  either  or  both  of  which  may  be  either 
circular  arcs  or  straight  lines,  in  principle.  The  most  interest¬ 
ing  applications  of  this  capability  probably  involve  the 
creation  of  "S-blades,"  composed  of  two  circular  arcs,  and 
"J-blades"  composed  of  a  straight-line  and  circular-arc 
combination . 

The  specification  of  this  two-segment  camber  line  requires 
the  boundary  condition  of  Equation  (11)  as  well  as  the  three 
boundary  conditions  necessary  for  the  circular-arc  camber  line. 

The  circular-arc  portions  are  computed  in  the  manner  of  the 
subsection  "Circular-Arc  Camber  Line"  using  the  inflection 
angle  and  the  inlet  (or  outlet)  angle.  Each  segment  is  scaled, 
and  positioned  with  respect  to  the  origin  (if  necessary).  A 
straight-line  segment  is  created  when  the  angle  at  the  inflection 
point  is  identical  to  the  inlet  (or  outlet)  angle. 

One  other  feature  of  the  programming  associated  with  the 
multiple-circular-arc  camber  line  should  be  noted:  an  ordinary 
circular-arc  camber  line  will  be  produced  if  the  inflection 
point  is  specified  at  the  blade  leading  edge.  Since  the  standard 
thickness  distribution  (described  in  Section  III)  is  applied 
to  the  multiple-circular-arc  camber  line,  this  feature  enables 
the  program  user  to  define  a  blade  using  the  circular-arc  camber 
line  and  either  the  standard  or  double  circular-arc  thickness 
distribution . 

Figure  11  shows  various  types  of  camber  lines  which  may  be 
specified  using  the  multiple-circular-arc  camber  line  option. 


Figure  12  shows  camber  lines  for  similar  "S-blade"  sections, 
illustrating  the  effect  of  varying  the  location  of  the 
inflection  point.  Camber  lines  for  three  "J-blade"  sections 
with  different  inflection  locations  are  shown  in  Figure  13. 


SECTION  III 


THE  SECTION  THICKNESS  DISTRIBUTION 


1.  STANDARD  THICKNESS  DISTRIBUTION 

For  the  same  reasons  that  it  was  desirable  to  define  the 
camber  line  in  equation  form,  it  is  also  practical  to  define 
the  thickness  distribution  this  way.  Five  features  were 
particularly  desired  for  a  supersonic  compressor  airfoil.  First, 
the  maximum  thickness  should  be  arbitrary.  Second,  it  should 
be  possible  to  locate  the  maximum  thickness  at  any  point  on 
the  downstream  half  of  the  airfoil.  Third,  there  should  be  no 
discontinuity  in  curvature,  meaning  that  derivatives  through 
the  second  should  be  continuous.  Fourth,  for  mechanical  reasons, 
the  airfoil  should  have  a  thickness  distribution  which  is 
continuously  convex;  i.e.,  the  second  derivative  should  not 
change  sign.  Fifth,  the  airfoil  should  have  the  slenderest 
leading  edge  that  is  compatible  with  the  first  four  conditions. 

Numerous  equations  were  examined  w ith  these  objectives  in 
mind.  However,  no  single  equation  was  found  which  seemed  to 
satisfy  all  of  the  foregoing  requirements  to  a  satisfactory 
extent.  The  method  finally  selected  consisted  of  defining 
the  thickness  distribution  with  two  equations:  one  from  the 
leading  edge  to  the  point  of  maximum  thickness  and  one  from 
there  to  the  trailing  edge.  At  the  point  of  juncture,  the 
thickness  and  the  first  and  second  derivatives  of  thickness 
were  equated.  In  order  to  prevent  a  reflex  curvature  from 
occurring  in  the  thickness  distribution  near  the  leading  edge, 
the  second  derivative  of  thickness  was  set  equal  to  zero  at  the 
leading  edge.  Thickness  was  considered  to  be  distributed 
normal  to  a  camber  line  of  unit  length.  Maximum  thickness  was 
located  at  a  point  Z  on  this  axis  which,  it  should  be  noted,  is 
not  the  chord  line.  The  thickness  of  the  leading  and  trailing 
edges  is  independently  specified  so  that  it  need  not  be  the 
same  at  both  edges.  At  the  leading  edge,  the  blade  surface  is 
completed  with  a  circular  arc.  At  the  trailing  edge,  the  blade 
was  truncated  by  connecting  the  two  end  points  with  a  straight 
line . 


The  boundary  conditions  chosen  for  the  two  segments  of 
the  thickness  distribution  may  be  summarized  as  follows. 
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For  the  Leading  Portion: 


At  x  -  0,  y  =  y( 
y"  =  0 


At  x  = 


z.  y  =  T/2 

y'  =  o 


For  the  Trailing  Portion: 


At  x  =  Z,  y  =  T/2 


(40) 


y'  =  0 


y "  =  y" 


(leading  portion) 


At  x  =  1 .0,  y  =  yt 


Note  that  y  represents  the  half-thickness  of  the  blade  element 
and  T  is  the  total  value  of  the  maximum  thickness.  Both  are 
scaled  in  proportion  to  a  camber  line  length  of  unity.  The 
parameter  Z  denotes  the  position  where  maximum  thickness  occurs 
on  the  camber  line.  For  the  leading  section,  a  simple  third- 
order  polynomial  was  used  in  the  form 

y  =  ax3  +  bx2  +  cx  +  d  (4l! 

y'  =  3ax2  +  2bx  +  c 
y"  =  6ax  +  2b 

Applying  the  boundary  conditions  for  this  portion  produces  the 
following  result  valid  from  0  to  Z: 


T/2  -  y, 


a  =  - 


3 ( T/ 2  -  y0) 


d  =  y. 
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Also,  at  x 

y" 


3(T/2  -  y0) 
_ 


(45) 


which  is  used  a?  a  boundary  condition  for  the  trailing  portion 
of  the  blade.  A  similar  third-order  equation  was  used  to 
define  the  trailing  portion  of  the  thickness  distribution.  This 
was  arranged 

y  =  e(x-Z) 3  +  f  ( x-Z) 2  +  g(x-Z)  +  h  (1|6) 

y*  =  3e(x-Z) 2  +  2f (x-Z)  +  g 
y"  =  6e(x-Z)  +  2f 


Applying  the  boundary  conditions  relating  to  the  trailing  portion 
leads  to  the  following  result  valid  from  x  =  z  to  1.0. 

3  ( T/2-y0 )  (T/2-y1)  (n) 

2Z2 ( 1 -Z)  (1-Z)3 

3(T/2-y0)  (48) 

f  = - - - 

2Z2 

9=0  (H9) 


h  =  T/2 


As  the  equations  for  the  two  portions  of  the  thickness 
distribution  are  cubics,  potentially  there  is  the  possibility 
of  introducing  inflection  points  into  the  thickness  distribution. 
One  boundary  condition  applied  to  the  equation  describing  the 
forward  portion  of  the  thickness  was  that  the  second  derivative 
be  zero  at  the  leading  edge.  Hence,  because  a  cubic  has  only 
one  point  of  inflection,  the  thickness  distribution  forward  of 
the  point  of  maximum  thickness  cannot  contain  an  inflection. 
Examination  of  the  equation  for  the  second  portion  of  the 
thickness  distribution  yields  the  result  that  there  is  a  minimum 
value  of  Z  (the  point  of  maximum  thickness)  such  that  the  in¬ 
flection  point  of  this  equation  is  not  on  the  blade  surface. 

This  is  given  by 
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1 


(50) 


Z  . 
mm 


1  + 


t 

y—0 


/T 

/T 


Thus  if  the  leading  and  trailing  edge  thicknesses  are  equal, 

Z  .  is  0.5,  and  will  increase  as  the  ratio  of  the  trailing 
eSge  to  leading  edge  thicknesses  increases. 

Figure  1*1  shows  thickness  distributions  for  three  values 
of  Z,  with  the  leading  and  trailing  edge  thicknesses  equal  in 
each  case. 

2.  DOUBLE-CIRCULAR-ARC  THICKNESS  DISTRIBUTION 


The  determination  of  the  double-circular-arc  thickness 
distribution  is  essentially  a  geometrical  exercise  involving 
the  translation  and  rotation  of  coordinate  systems  and  the  use 
of  polar  coordinates  for  the  sake  of  simplicity  and  convenience. 


The  most  desirable  coordinate  system  would  define  the 
circular  arc  camber  line  as  part  of  a  circle  centered  at  the 
origin,  such  that  the  camber  line  represents  the  portion  of  the 
circle  in  the  range 


TT 

2 


-  <  ill  < 

2  “ 


!L  +  1 
2  2 


(51) 


when  i ji  is  measured  counterclockwise  from  the  positive  x  axis 
and,  from  geometry, 


*  =  ai"a2 


(52) 


This  combined  translation-rotation  is  accomplished  by  the 
coordinate  change 

x'  =  ( x  -  x  j )  cos  p  +  (y-yj  sin  p  (53) 

y'  =  -  ( x- x j )  sin  p  +  (y-yj  cos  p 


where 


(5*0 
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From  the  geometry  of  the  camber  line. 


W  =  sin-1 


<x0  -  D 


In  the  new  coordinate  system,  then, 


x'2  +  y'2  =  R2 


In  this  system,  either  surface  of  the  airfoil  may  be  expressed 
as  a  portion  of  a  circle  centered  on  the  y*  axis.  That  is, 

x’2  +  (y1  -  y  )2  =  R  2 

J  S  s 


or,  expanding, 


x'!  +  y ' 2  -  2y '  ys  +  y$2  =  Rs2 


In  cylindrical  polar  coordinates.  Equations  (56)  and  (58) 
may  be  expressed  as  follows 


r!  ■  R 

r22  +  ys2  -  2  r2  ys  sin  f  -  Rs2 


where  the  relevant  range  of  ip  is  as  given  in  Equation  (51)- 

The  desired  thickness  is  the  difference  (Ar)  between  the 
expressions  for  the  two  arcs.  First,  however,  Equation  (60) 
must  be  solved  for  r2  rather  than  r2  .  Using  the  binomial 
theorem, 


r2  =  ys  sin  ip  ±  J  Rg2  -  yg2  cos2  ip 


If  it  is  stipulated  that  Equation  (6l)  represents  the  surface 
of  the  blade  which  lies  further  from  the  origin  than  does  the 
camber  line,  then  the  thickness  t  in  a  positive  sense  is  given 
by  the  expression  r2  -  rt . 
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t 


y  s  i  n  ^  ± 

3 


C0S2iJ, 


R 


(62) 


The  proper  sign  of  the  second  term  (+)  is  chosen  by 
requiring  that  at  =  ir/2, 

t  =  y  +  R  -  R  (63) 

s  s 


Therefore,  Equation  (62)  becomes 


t  =  y  sin  ip  +  J  R  2  -  y  2  cos2if>  -  R 

S  v  S  S 


(6H) 


A  further  rotation  of  the  coordinate  system  such  that 


(65) 


yields 


t  =  y  cos 


[M  +V 


R  2  -  y  2  si 

s  s 


•  2  r  4> 

in  |  -  v 


-n 


(66) 


-  R 


in  the  range  OS?  u  ^  <j>. 

Expressing  Equation  (66)  in  terms  of  the  normalized  arc 
length,  s. 


t  =  y  cos 
J  s 


2  -  s  )  * 


V 


+  VRS2  -  ys2  sin2 


(67) 


2  ‘  s  )  * 


-  R 


in  the  range  0  <  s  <  1. 

The  expressions  for  ys  and  Rs  in  Equation  (67)  can  be 
determined  geometrically  in  terms  of  T/2  and  rQ,  the  edge  radius 


R  =  rA  + 
s  0 


R  sin  4> 


(68) 


sin 


Arccos 


l-A' 
1+ A2 
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yc  =  R  -  R  +  i- 

s  s  2 


(69) 


where 


(T/2  -  rQ)  +  R  ( 1  -  cos  <f>/2) 
R  sin  <p/2 


(70) 


Equations  (68),  (69)  and  (70)  yield  all  the  unknown 
constants  of  the  thickness  distribution  in  terms  of  known 
quantities. 
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SECTION  IV 


CARTESIAN  COORDINATES  FOR  THE  BLADE 


The  two  previous  sections  have  described  the  methods 
derived  to  specify  individual  blade  sections.  When  located  as 
desired  relative  to  the  blade  stacking  axis,  the  section 
coordinates  are  the  coordinates  of  the  streamsurface  blade 
section.  A  series  of  such  sections  on  all  given  stream  surfaces 
specifies  the  envelope  of  the  blade,  but  the  surface  coordinates 
are  not  in  a  convenient  form  for  manufacturing  purposes.  A 
method  of  determining  the  coordinates  of  plane  sections  through 
the  blade  at  a  series  of  locations  along  the  stacking  axis  is 
now  described. 


The  method  relies  upon  the  use  of  the  "spline-curve"  for 
interpolation  (or  extrapolation)  of  the  coordinates  of  the 
blade  surfaces  and  therefore  the  theory  of  the  curve  is  given. 

The  same  interpolation  scheme  is  also  used  elsewhere  in  the 
computer  program. 

Given  the  Cartesian  coordinates  x  and  y  for  a  series  of 
points,  the  problem  is  to  find  a  curve  y  =  f(x)  that  passes 
through  each  given  point.  A  series  of  algebraic  cubics  is 
used,  one  equation  applying  between  each  adjacent  pair  of 
points.  The  four  coefficients  for  each  cubic  permit  the 
matching  of  the  first  and  second  derivatives  of  the  two  equations 
applying  at  each  point  (except  the  first  and  last),  in  addition 
to  passing  the  curves  through  the  specified  points.  Two 
additional  cu  ditions  are  required  to  establish  the  coefficients. 
Those  chosen  here  are  that  the  second  derivatives  of  the 
equations  at  the  first  and  last  points  are  zero.  Other 
conditions  are  possible,  but  less  appropriate  for  the  current 
purpose.  A  curve  of  the  form  described  is  that  theoretically 
assumed  by  a  thin,  uniforinly-flexible  beam  pinned  at  each 
specified  point,  hence  the  name  "spline-curve." 


By  applying  the  conditions  described  above,  the  following 
system  of  equations  may  be  produced.  In  the  interval  between 
any  two  adjacent  points,  such  that  xn_,<x<xn,  the  equation 
of  the  spline-curve  is  !i 


y  = 


1  .  (x  -x)3  M  (x-x  )3 

n~1  n _  +  —  - Uli —  +  (x  -x) 

6  <Wi>  6  <Wi> 


r  ,  ii  . 

n  - 1  n  - 1 


X-X  , 
n  n- 1 


(x-x  n) 
'  n-1 


X-X  . 
n  n-1 


—  (x  -X  - ) 
e.  '  n  n-1 


<VXn-l> 


(71) 
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'ind  the  s 3 ope  of  the  spline-curve  In 


r 


dy 

dx 


M 


(x-x)2  M  (x-x  .)2  /y  -y 
n-1  n  +  _H  n-1'  (Jn  •'n-l 


(72) 


M  -M  . 
n  n-1 


2  ( x- x  )  2  (x-x  .)  \x-x  . 

n  n-17  '  n  n-1'  ' • n  n-1 


)  <Wl> 


The  coefficients  M  are  given  by 

total  number  of  specified  points,  and  N 


=0.0  where  N  is  the 


M 


0  -B  M  . 
n  n  n  +  1 


(73) 


A:  =  1  .0 


B  =  0.0 


D1  =  0.0 


(x  -x  ,)  (x-x  ,)B 
_  v  n+1  n-1'  n  n-1  n-1 


n-1 


X  -X 
n+1  n 


yn.ryn  Wl  ([ V V d  Dn- : 


X  ^  -X  X-X  . 

n+1  n  n  n-1 


6  A 


n-1 


Hence  the  coefficients  M  may  be  found  by  first  computing  the 
coefficients  A,  B  and  D. 

In  order  to  use  the  spline-curve  for  interpolation,  it  is 
convenient  first  to  compute  the  M  coefficients  for  all  the 
specified  points,  and  then  to  evaluate  the  value  of  the  function 
y  at  the  desired  value  (or  values)  of  x,  using  the  coefficients 
appropriate  to  the  interval  in  which  x  lies.  Should  extrapolation 
outside  the  specified  range  of  values  of  x  be  required,  this  is 
-lone  assuming  a  straight  line  at  the  slope  of  the  spline-curve 
at  the  appropriate  end-point. 
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The  determination  of  the  Cartesian  coordinates  of  the 
manufacturing  sections  involves  two  steps.  First,  Cartesian 
coordinates  for  the  points  defining  the  streamsurfaces  are 
obtained.  Second,  coordinates  for  the  manufacturing  sections 
are  interpolated  (or  extrapolated)  at  the  desired  plane  sections 
through  the  blade. 

Figure  15  illustrates  the  geometry  involved  in  the  deter¬ 
mination  of  the  Cartesian  coordinates  of  the  streamsurface 
sections.  The  radius  of  the  streamsurface  is  defined  at  a 
number  of  axial  locations.  A  spline-curve  is  passed  through 
these  points,  and  the  slope  of  the  streamsurface  is  calculated 
(as  the  slope  of  the  spline-curve)  at  100  points  distributed 
uniformly  on  the  x-axis  between  the  first  and  last  given  points. 
Then  an  m-x  table  is  constructed  by  numerically  integrating  for 
the  streamsurface  using 

(7*0 


The  m  coordinate  of  any  point  on  the  streamsurface  blade  section 
is  the  "axial"  coordinate  of  the  point  derived  from  the  camber 
line  and  thickness  distribution.  Thus  linear  interpolation  of 
the  m-x  table  yields  the  corresponding  x  coordinate.  Use  of 
a  100-point  table  assures  that  linear  interpolation  is  adequate. 
The  r-x  table  given  enables  the  radius  of  the  streamsurface  at 
any  value  of  x  to  be  found.  Spline-curve  interpolation  is  used 
for  this  purpose. 

The  "y"  coordinate  of  any  point  derived  from  the  camber 
line  and  thickness  distribution  must  now  be  related  to  the 
circumferential  direction  on  a  streamsurface.  This  can  be 
accomplished  by  determining  the  angle  e  (defined  in  Figure  15) 
in  the  following  manner.  First,  the  origin  of  the  m-system  is 
shifted  so  that  the  origin  (m  =  0)  corresponds  to  the  location 
of  the  centroid  of  the  section  in  the  cascade  plane.  The 
camber  line  slopes  in  the  cascade  plane  must  then  be  reproduced 
on  a  streamsurface.  This  requires  that  at  each  point 


de  _  d,y 
dm  ”  dx 


y 1 


(75) 


It  is  possible  to  determine  the  variation  of  e  with  m  by 
integrating  Equation  (75)  with  respect  to  m. 


2*) 


(76) 


where  both  the  streamsurface  radius  and  the  camber  line  slope 
are  known  functions  of  x,  hence  m. 

At  m  =  0, 


e 


o 


(77) 


where  is  the  distance  between  the  section  centroid  and  the 
camber  line  measured  in  the  y  direction  in  the  cascade  plane, 
and  r^  is  the  streamsurface  radius. 

Therefore,  from  Equations  (7 6)  and  (77), 


e(m)  = 


(78) 


For  each  point  on  the  blade  surface  in  the  cascade  plane, 
the  appropriate  e  is  determined  through  a  six  step  procedure: 

(1)  m  is  determined  from  the  m-x  table 


(2)  G camber  line  is  determined  from  Equation  (78) 

<3LAiiS„?“ned  fr°m  y  =  ^surface  '  Camber  line 
m  the  cascade  plane 

(*0  r  is  determined  from  the  value  of  m  in  conjunction 
with  the  m-x  and  r-x  tables 


(5)  Ae  is  determined  from  As  =  Ay/r 
^  Gsurface  -  ecamber  line  +  ^G 

The  two  remaining  Cartesian  coordinates  of  the  point  then 
follow  from 


z  =  r  cos  e 
and 

y  =  r  sin  e 


(79) 


The  interpolation  of  x  and  y  coordinates  for  the  manufactur¬ 
ing  sections  (on  constant-z  planes)  is  achieved  by  fitting 
the  spline-curve  through  corresponding  points  on  each  stream- 
surface  section.  The  same  number  of  points  is  used  to  define 
each  blade  section.  Hence,  the  spline-curve  is  used  first 
to  fit  x  as  a  function  of  z  for  the  first  point  on  the  suction 
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surface  of  each  section.  The  value  of  x  is  then  interpolated  at 
each  desired  value  of  z.  Then  this  is  repeated  for  the 
coordinate  y,  and  then  for  all  other  corresponding  points 
describing  the  streamsurface  sections.  Thus  an  equal  number  of 
similarly  distributed  points  describing  the  manufacturing 
sections  is  obtained. 


SECTION  V 


SECTION  PROPERTIES 


The  stacking  axis  of  the  blade  is  passed  through  each 
streamsurface  section  either  at  its  leading  or  trailing  edge, 
or  at  a  point  specified  relative  to  the  centroid  of  the  section. 
Because  the  streamsurface  sections  are  in  general  not  planar, 
the  centroids  of  the  manufacturing  sections  will  not  generally 
lie  precisely  on  the  stacking  axis  when  the  streamsurface 
sections  are  stacked  about  their  centroids.  By  determining  the 
locations  of  the  centroids  of  the  manufacturing  sections  so 
obtained,  it  is  possible  to  estimate  the  offsets  that  must 
be  applied  when  restacking  the  streamsurface  sections  in  order 
that  the  centroids  of  the  manufacturing  sections  will  be  located 
as  desired  relative  to  the  stacking  axis.  Thus  there  is  a 
requirement  to  determine  the  locations  of  the  centroids  of 
both  the  streamsurface  and  manufacturing  sections  if  the  blade 
is  to  be  stacked  on  the  centroids. 

To  further  assist  in  any  mechanical  analysis  of  the  blade, 
the  areas,  second  moments  of  area,  principal  axes,  and 
principal  second  moments  of  area  for  both  the  streamsurface 
and  manufacturing  sections  are  also  determined.  The  methods 
used  to  determine  these  quantities  are  given  below. 

1.  STREAMSURFACE  SECTIONS 

The  location  of  the  centroid  of  an  area  relative  to  an 
arbitrary  origin  is  given  by 


In  the  computer  program  these  integrals  are  performed 
numerically  by  considering  the  section  to  be  composed  of  a 
semicircle  formed  by  the  leading  edge  and  a  series  of 
quadrilaterals.  The  quadri laterals  are  defined  by  the  series 
of  points  defining  the  camber  line  and  the  values  of  the 
section  thickness  at  these  points,  as  shown  in  Figure  16.  The 
origin  for  the  calculation  is  the  leading  edge  of  the  blade  and 
each  quadrilateral  is  assumed  to  be  a  rectangle.  The  length 
of  one  side  is  the  camber  line  length  enclosed  by  the 
quadrilateral,  and  the  length  of  the  other  is  the  mean  section 
thickness  at  the  ends  of  the  quadrilateral.  Hence  the  above 
relationships  are  evaluated  as 


?7 


rxdA  = 


(81) 

!n±l  (s  _s  )  (tn<-i't") 

' sn+l  V  2 


K  Od) 

f„  ^o2  4*»  .  Y'  /ymn+ymn+l  .  ,  (V‘n+11) 

yydA  ■  ■  2  3 it  1  1'^1  \  2  *sn+rsn^  2  j 


MA  = 
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For  the  case  of  the  double-circular-arc  blade,  a  term 
representing  the  contribution  of  the  trailing  edge  semicircle, 
similar  to  the  term  for  the  leading  edge  contribution,  is  added 
to  each  of  Equations  (81),  (82)  and  (83). 

The  second  moments  of  area  I  and  I  ,  and  products  of 
inertia,  I  ,  are  found  by  again  considering  the  section  to  be 
composed  ofxa  series  of  rectangles,  as  shown  in  Figure  16. 

(The  contribution  of  the  edge(s)  is  neglected.)  The  principal 
second  moments  of  one  rectange  about  its  centroid  are 


1 X ‘  =  ^sn+l~sn^ 


tn+l+tn 


■  (^)  <vrsn>3  /’2 


The  mean  camber  line  angle,  C,  at  the  two  ends  of  the 
rectangle  is  used  to  determine  the  second  moments  of  area  of 
the  rectangle  about  its  centroid  from 

^  y  1  -  ^  v  '  ( 

Ix  -  -x  2  y  +  -x  i  y  cos (2C) 


IY.+IU.  iv.-i  . 

■  -X  ~2  X- - "X  2  ^  •  CO s  ( 2 C ) 
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The  product  of  inertia  of  the  rectangle  is  also  determined  from 


I 


xy 


(Ix-Iy)  cos ( 2C)  -  Ix,+Iy, 
2  s  i  n  ( 2  C ) 


(88) 


The  contributions  of  the  rectangle  to  the  second  moments  of 
area  of  the  section  about  its  centroid  are  given  by 


^2  (89) 

Al  (section  centroid)  =  I  +  A(y-y) 

a  X 


2  (90) 

Al  (section  centroid)  =  I  +  A(x-x) 

J  j 


where  A  is  the  area  of  the  rectangle,  and  x  and  y  are  the 
coordinates  of  the  center  (that  is,  centroid)  of  the  rectangle. 
Also,  the  contribution  of  the  rectangle  to  the  product  of 
inertia  of  the  section  about  its  centroid  is  given  by 

Al  (section  centroid)  =  I  +  A(x-x)  (y-y)  (91) 

xy  xy 


The  summation  of  all  such  contributions  yields  the  second 
moments  of  area  and  product  of  inertia  of  the  section  about  its 
centroid. 

The  inclination  of  the  principal  axes  of  the  section,  B, 
is  given  by 


B 


2  arctan 


(92) 


The  principal  second  moments  of  area  of  the  section  about  its 
centroid  are  given  by 


I  +1  I  -I 

i.._  =  —  cos ( 2B )  -  I  s i n ( 2 B ) 


(93) 
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xy 


i+i  I  -i 

Iyp  =  -V^  -  “V*  C0S(2B)  +  !xy  s i n ( 2 B ) 
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The  various  calculations  described  above  are  made  for  a 
"normalized"  streams urface  blade  section  (that  is  one  having  an 
axial  chord  of  unity),  and  are  subsequently  scaled  as  appropriate, 
according  to  the  specified  chord  length. 

2.  MANUFACTURING  SECTIONS 

The  methods  used  to  determine  the  properties  of  the 
manufacturing  sections  are  essentially  the  same  as  those  described 
above  for  the  streamsurface  sections.  However,  in  this  case 
the  section  geometry  is  described  by  the  coordinates  of  a  series 
of  points  on  the  two  blade  surfaces  as  shown  in  Figure  16. 

The  section  is  again  considered  to  be  composed  of  a 
semicircle  at  the  leading  edge  and  a  number  of  rectangles. 

The  radius  of  the  leading  edge,  yQ,  is  calculated  to  be  a  half 
of  the  distance  between  the  first  points  on  the  blade  pressure 
and  suction  surfaces.  The  angle  of  the  blade  section  at  the 
leading  edge  is  obtained  from  the  first  two  points  on  each 
blade  surface,  thus 


/yp2+ys2-yp1-ys1  \  (95) 

(vxp2+xs2-xp1-xs1  J 


The  lengths  of  the  sides  of  each  assumed  rectange  are  determined 
from  the  coordinates  of  the  points  at  the  four  corners,  thus 
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Hence,  the  section  area  and  the  location  of  the  centroid  are 
determined  from 
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As  for  the  st.  reams  ur  face  sections,  an  aiditional  term  is 
inserted  in  each  of  Equations  (98),  (99)  and  (100)  for  the 
double-circular-arc  blade,  to  account  for  the  contribution  of 
the  trailing  edge  semi-circle. 


The  second  moments  of  area  and  product  of  inertia  of  the 
section  are  again  found  by  first  determining  the  principal 
second  moments  of  each  rectangle  about  its  centroid,  and  then 
determining  the  contribution  of  each  rectangle  to  the  second 
moments  and  product  of  inertia  of  the  section.  Referring  again 
to  Figure  16 ,  we  have  for  one  rectangle 


Ix,  =  ml  1 12 


(101) 


Iyl  =  Zm 3 / 1 2 


(102) 


In  order  to  determine  the  second  moments  of  area  of  the 
rectangle  about  is  centroid  In  directions  parallel  to  the  axis 
of  the  section  coordinates,  the  angle  of  inclination  of  the 
rectangle  is  found  from 
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arctan 
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Then  the  second  moments  of  area  and  product  of  inertia  for 
the  rectangle  are  found  using  Equations  (86),  (8?)  and  (88) 
arain,  and  the  contributions  of  the  rectangle  to  the  properties 
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where 
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and  du,  t,  and  S  are  the  elemental  lancr-t-v-i  ~-i  <-u  ■  . 
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SECTION  VI 


BLADE  CHARACTERISTICS 


1.  BLADE  VOLUME  CALCULATION 

Except  when  IPRINT  =  1,  a  calculation  of  the  volume 
enclosed  by  the  blade  between  the  innermost  and  outermost 
streamsurfaces  is  made.  This  is  done  from  the  Cartesian 
coordinates  of  the  points  describing  the  streamsurface  blade 
sections.  The  total  volume  is  obtained  by  summing  the  volumes 
enclosed  between  each  adjacent  pair  of  streamsurfaces.  These 
volumes  in  turn  are  obtained  by  summing  two  types  of  sub¬ 
volume.  First,  there  is  the  volume  enclosed  by  the  leading 
edge  radius.  This  is  determined  by  finding  the  mean  area  of 
the  two  leading  edge  semicircles  on  the  two  streamsurfaces, 
and  multiplying  by  the  mean  height  at  the  two  ends  of  the  semi¬ 
circle  diameters.  The  diameter  of  the  leading  edge  is  calculated 
to  be  the  distance  between  the  first  points  on  the  section 
pressure  and  suction  surfaces  using  the  streamsurface  Cartesian 
coordinates.  Thus  the  areas  computed  are  those  projected  onto 
a  constant-z  plane,  and  the  volume  is  given  directly  by 
multiplying  by  the  mean  z  dimension.  Second,  there  is  a 
series  of  volumes  enclosed  by  corresponding  quadrilaterals  on 
the  two  streamsurfaces.  Again,  the  sub -volume  is  determined 
from  the  mean  of  the  areas  of  the  two  quadrilaterals  projected 
onto  a  constant-z  plane,  and  the  mean  distance  between  them. 

The  quadrilaterals  are  each  defined  by  two  adjacent  points 
on  the  section  suction  surface  and  the  two  corresponding 
points  on  the  pressure  surface.  The  area  of  the  quadrilaterals 
is  obtained  by  assuming  them  to  be  rectangles  of  side  Z  and  m, 
as  defined  in  Equations  (96)  and  (97),  except  that  in  this 
case  Cartesian  coordinates  of  the  streamsurface  sections  are 
used.  The  distance  between  the  two  rectangles  is  taken  to  be 
the  mean  of  the  distance  between  the  four  corresponding  corners 
of  the  rectangle.  Thus  the  total  volume,  V,  is  given  by 
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where  J  is  the  total  number  of  streamsurfaces,  z  is  the  mean 
coordinate  at  each  end  of  the  leading  edge  semicircle  diameter, 

A  is  the  mean  of  the  two  leading  edge  semicircle  areas,  N  is 
the  number  of  points  defining  each  blade  surface,  and  z  is  the 
mean  z-coordinate  at  the  four  corners  of  the  rectangle. n  Equation 
(106)  is  modified  to  contain  a  term  similar  to  the  first  to 
represent  the  trailing  edge  contribution  for  the  double-circular- 
arc  blade. 


33 


2. 


CALCULATIONS  FOR  AERODYNAMIC  ANALYSIS 


Certain  aerodynamic  analyses  of  the  blade  require  a 
description  of  the  blade  on  cylindrical  surfaces.  As  an  option, 
the  computer  program  will  produce  quantities  necessary  for  this 
description  at  a  specified  streamsurface-computing  station 
intersection. 

First,  the  angular  position  of  the  camber  line  with  respect 
to  the  stack  axis  at  a  streamsurface-computing  station  inter¬ 
section  is  specified  in  terms  of  $,  which  is  the  value  of 
ecamber  line  on  s breams ur face  at  that  particular  intersection. 

The  physical  passage  blockage  (b)  due  to  the  presence  of 
the  blades  is  determined  as  a  percentage  of  the  passage 
circumference  in  terms  of  the  number  of  blades  in  the  blade 
row,  N,  and  a,  the  angle  subtended  on  the  cylindrical  surface 
by  each  blade. 


b 


(107) 


The  blade  lean  angle,  £,  with  respect  to  the  radial 
direction  at  a  given  point  is  obtained  from  the  si  ope  of  a 
spline- curve  fit  through  the  y-z  Cartesian  coordinates  of  the 
streamsurface  section  camber  lines  at  the  particular  axial 
location.  Thus 

£  =  <t>  -  Arctan^0  (108) 


Two  other  quantities  are  needed  to  produce  the  proper 
mean-camber  line  angle  on  the  cylindrical  surface:  the  local 
computing  station  inclination,  y,  obtained  from  the  specified 
station  description;  and  the  local  streamsurface  inclination,  y, 
obtained  from  the  specified  r-x  streamsurface  description.  f 
Together  with  the  camber  line  angle  on  the  streamsurface,  a°: 


tan  a£ 


tan  y  tan  £  +  tan  a*/cos  y 
1-tan  y  tan  y 


(109) 


Equation  (109) 


is  taken  from  Reference  3. 
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SECTION  VII 


USE  OF  THE  COMPUTER  PROGRAM 


Basic  information  required  by  the  user  to  run  the  computer 
program  is  given  in  this  section.  The  various  input  data  items 
are  defined  first,  and  the  input  data  format  is  then  specified. 

A  description  of  the  output  data  that  may  be  expected  is  given. 
(Implementation  of  the  program  on  a  computing  system  is  not 
discussed  here,  but  in  the  section  entitled  "Computer  Program 
Details". ) 

1.  DEFINITION  OF  INPUT  DATA  ITEMS 

TITLE  An  alphanumeric  title  of  72  characters  that  may 

be  used  to  identify  a  run. 

NLINES  The  number  of  streamsurfaces  which  are  defined 

on  and  on  which  blade  sections  will  be  designed. 
Must  satisfy  2  <  NLINES  <15- 

NSTNS  The  number  of  computing  stations  at  which  the 

streamsurface  radii  are  specified.  Must  satisfy 
3  <  NSTNS  <  10  . 

NZ  The  number  of  constant-z  planes  on  which 

manufacturing  (Cartesian)  coordinates  for  the 
blsde  are  required.  Must  satisfy  3<NZ<15. 

NSPEC  Tne  number  of  radially-disposed  points  at  which 

the  parameters  of  the  blade  sections  are 
specified.  Must  satisfy  1<  NSPEC  <15. 

NPOINT  The  number  of  points  that  will  be  generated  to 

specify  the  pressure  and  suction  surfaces  of 
each  blade  section.  Must  satisfy  2<NPOINT<80. 
Generally,  no  less  than  30  should  be  used.  It 
will  be  advantageous  to  specify  80  points  when 
precision  plots  of  the  sections  are  to  be 
produced  directly  by  the  program. 

NBLADE  The  number  of  blades  in  the  blade  row. 

ISTAK  If  ISTAK  =  0,  the  blade  will  be  stacked  at  the 

leading  edge . 

If  ISTAK  =  1,  the  blade  will  be  stacked  at  the 
trailing  edge  . 

If  ISTAK  =  2,  the  blade  will  be  stacked  at,  or 
offset  from,  the  section  centroid. 
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I PUNCH 


ISECN 


IFCORD 


IFPLOT 


An  integer  which  indicates  whether  the  quantities 
necessary  for  aerodynamic  analysis  of  the  blade 
are  to  be  produced  on  punched  cards . 

If  ISECN  =  0,  the  blade  will  be  constructed 
using  the  polynomial  camber  line  and  the 
standard  thickness  distribution. 

If  ISECN  =  1,  the  exponential  camber  line  and 
the  standard  thickness  distribution  will  be 
used. 

If  ISECN  =  2,  a  blade  using  the  circular  arc 
camber  line  and  the  double-circular-arc 
thickness  distribution  will  be  used. 

If  ISECN  =  3j  the  standard  thickness 
distribution  and  any  of  the  camber  lines 
described  in  the  sub-section,  "Multiple 
Circular-Arc  Camber  Lines,"  may  be  used. 

If  IFCORD  =  0,  the  meridional  projections  of 
the  streamsurface  blade  section  chords  are 
specified. 

If  IFCORD  =  1,  the  streamsurface  blade  section 
chords  are  specified. 

Where  CALCOMP  software  is  incorporated  into 
the  computing  system,  IFPLOT  specifies  the 
creation  of  precision  plots.  (Further  information 
regarding  the  requirements  for  this  are  given 
in  the  section  entitled  "Computer  Program 
Details". ) 

If  IFPLOT  =  0,  no  plots  will  be  produced. 

If  IFPLOT  =  1,  a  plot  of  the  streamsurface 
sections  will  be  produced.  All  NLINES  sections 
are  shown  superimposed.  The  origin  for  each 
section  plot  is  offset  from  the  centroid  of 
the  section  by  distances  specified  by  DELX 
and  DELY. 

If  IFPLOT  =  2,  a  plot  of  the  manufacturing 
sections  will  be  produced.  The  origin  is  the 
blade  stacking,  axis,  and  all  NZ  sections  are 
shown  superimposed. 

If  IFPLOT  =  3,  both  of  the  plots  described  for 
IFPLOT  =  1  and  2  will  be  produced. 
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if  JPP1.0T  =  -I,  individual  plots  of  each  of  the 
manufacturing  sections  will  be  produced.  The 
axes  are  rotated  clockwise  by  the  section 
staffer  angle  for  each  plot. 

The  input  data  is  always  listed  by  the  program. 
Details  of  the  s t reams. ur face  and- manufacturing 
suctions  are  printed  as  prescribed  by  IPRINT. 

if  I  PRINT  =  0,  details  of  streamsurface  and 
manufacturing  sections  are  printed. 

if  IPRINT  =  1,  details  of  streamsurface 
sections  are  printed. 

If  IPRINT  =  2,  details  of  the  manufacturing 
sections  are  printed. 

The  NZ  manufacturing  sections  are  equispaced 
between  z  equals  "INNER  and  ZOUTER. 

When  precision  plots  are  produced,  SCALE  is 
the  scale  factor  employed.  , 

Tills  is  the  axial  coordinate  of  the  stacking 
axis  for  the  blade,  relative  to  the  same 
origin  as  used  for  the  station  locations,  XSTA. 

The  size  (inches)  of  the  plotter  to  be  used 
in  the  creation  of  precision  plots. 

The  number  of  points  provided  to  specify  the 
shape  of  a  computing  station. 

If  KPTS  =  1,  the  computing  station  is  upright 
and  linear. 

If  KPTS  =  2,  the  computing  station  is  linear 
and  either  upright  or  inclined. 

If  KPTS  >  2,  a  spline  curve  is  fit  through  the 
points  provided  to  specify  the  shape  of  the 
station . 

If  IFANGS  =0,  the  calculations  of  the 
quantities  required  for  aerodynamic  analysis 
will  be  omitted  at  a  particular  computing 
station. 

If  IFANGS  =  1,  these  calculations  will  be 
performed  at  that  station. 
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XSTA  An  array  of  KPTS  axial  coordinates  (relative  to 

an  arbitrary  origin)  which,  together  with  RSTA, 
specify  the  shape  of  a  particular  computing 
station . 

RSTA  An  array  of  KPTS  radii  which,  together  with 

XSTA,  specify  the  shape  of  a  particular 
computing  station. 

R  The  streamsurface  radii  at  NLINES  locations 

at  each  of  the  NSTNS  stations. 

ZR  The  variation  of  properties  of  the  streamsurface 

blade  sections  is  specified  as  a  function  of 
streamsurface  number.  The  various  quantities 
are  then  interpolated  (or  extrapolated)  at 
each  streamsurface.  The  streamsurfaces  are 
numbered  consecutively  from  the  innermost 
outward,  starting  with  1.0.  ZR  must  increase 
monotonically ,  there  being  NSPEC  values  in  all. 

B1  The  blade  inlet  angle.  (Figure  17  illustrates 

the  geometry  of  a  typical  blade  section.) 

B2  The  blade  outlet  angle. 

PP  If  ISECN  =  0,  PP  is  the  ratio  of  the  second 

derivative  of  the  camber  line  at  the  leading 
edge  to  its  maximum  value.  Must  satisfy 
-2.0  <  PP  <  1.0. 

If  ISECN  =  1,  PP  is  the  ratio  of  the  second 

derivative  of  the  camber  line  at  the  leading 

edge  to  its  maximum  value  forward  of  the 
inflection  point.  Must  satisfy  0.0  <  PP  <  1.0. 

If  ISECN  =  2  or  3j  PP  is  superfluous. 

QQ  If  ISECN  =  0,  QQ  is  the  ratio  of  the  second 

derivative  of  the  camber  line  at  the  trailing 
edge  to  its  maximum  value.  Must  satisfy 
0.0  <  QQ  <  1.0. 

If  ISECN  =  1,  QQ  is  the  ratio  of  the  second 

derivative  of  the  camber  line  at  tl  e  trailing 

edge  to  its  maximum  value  rearward  of  the 
inflection  point.  Must  satisfy  0.0  <  QQ  <  1.0. 

MOTE:  Some  comments  regarding  useful  values  for 
PP,  QQ,  S  and  BS  are  made  in  the  section 
entitled  "Examples  of  Use  of  Computer  Program." 
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If  ISECN  =  2  or  3,  QQ  is  superfluous. 

The  ratio  of  blade  leading  edge  radius  to 
chord. 

TC  The  ratio  of  blade  maximum  thickness  to  chord. 

TE  The  ratio  of  blade  trailing  edge  half-thickness 

to  chord. 

If  ISECN  =  2,  TE  is  superfluous. 

Z  The  location  of  the  blade  maximum  thickness, 

as  a  fraction  of  camber  line  length. 

If  ISECN  =  2,  Z  is  superfluous. 

CORD  If  IPCORD  =  0,  CORD  is  the  meridional  projection 

of  the  blade  chord. 

If  IPCORD  =  1,  CORD  is  the  blade  chord. 

DELX,  The  stacking  axis  passes  through  the  stream- 

DELY  surface  blade  sections,  offset  from  the 

centroids,  leading,  or  trailing  edge  by  DELX 
and  DELY  in  the  x  and  y  directions  respectively. 

S,  BS  If  ISECN  =  1  or  3,  S  and  BS  are  used  to  specify 

the  location  of  the  inflection  point  (as  a 
fraction  of  the  meridionally-projected  chord 
length)  and  the  change  in  camber  angle  from 
the  leading  edge  to  the  inflection  point. 

If  the  absolute  value  of  the  angle  at  the 
inflection  point  is  larger  than  the  absolute 
value  of  Bl,  BS  should  have  the  same  sign  as 
Bl;  otherwise,  Bl  and  BS  should  be  of  opposite 
sign . 

2.  INPUT  DATA  FORMAT 

Data  input  is  by  punched  card,  and  three  formats  are  used. 
The  first  card  only  is  alphanumeric,  using  the  first  72  columns 
on  the  card.  This  is  followed  by  one  card  of  integers.  Integers 
are  each  placed  in  a  field  of  three  locations,  starting  in 
Column  1  of  the  card.  No  decimal  point  may  be  used,  and  the 
number  must  occupy  the  right-most  locations  of  the  allocated 
field.  Real  numbers  are  each  placed  in  fields  of  12  locations, 
starting  in  Column  1  of  the  card.  Decimal  points  should  be 
included  to  ensure  correct  interpretation  of  the  data,  and  the 
numbers  may  be  placed  anywhere  in  the  allocated  fields. 
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In  the  following  chart,  one  line  corresponds  to  one  card. 


TITLE 

NLINES  NSTNS  NZ  NSPEC  NPOINT  MBLADE  ISTAK  IPUNCH  ISECN  IFCORD  IFPLOT  IPRINT 


ZINNER  ZOUTER  SCALE  STACKX  PLTSZE 
KPTS  I FANGS 

XSTA  RSTA  }  repeated  KPTS  times  ) 

?  repeated  NSTNS  times 
R  f  repeated  NLINES  times  ) 

ZR  B1  B2  PP  QQ  RLE  1 

TC  TE  Z  CORD  DELX  DELY  [  repeated  NSPEC  times 

S  BS  -  if  ISECN  =  1  or  3  only  / 

Listings  of  sample  input  data  decks  are  included  under 
"Examples  of  Use  of  Computer  Program." 

3.  OUTPUT  DATA 


Printed  output  from  the  program  may  be  considered  to 
consist  of  four  sections:  a  printout  of  the  input  data,  details 
of  the  blade  sections  on  each  streamsurface ,  a  listing  of 
quantities  required  for  aerodynamic  analysis,  and  details  of  the 
manufacturing  sections  determined  on  the  constant-z  planes.  These 
are  briefly  described  below.  In  the  explanation  which  follows, 
parenthetical  statements  are  understood  to  refer  to  the  particular 
case  of  the  double-circular-arc  blade  (ISECN  =2). 

The  input  data  printout  includes  all  quantities  read  in, 
and  is  self-explanatory. 

Details  of  the  streamsurface  blade  sections  are  printed  if 
IPRINT  =  0  or  1.  Listed  first  are  the  parameters  defining  the 
blade  section.  These  are  interpolated  at  the  streamsurface  from 
the  tables  read  in.  Then  follow  details  of  the  blade  section  in 
"normalized"  form.  The  blade  section  geometry  is  given  for  the 
section  specified,  except  that  the  meridional  project  of  the 
chord  is  unity.  For  this  section  of  the  output,  the  coordinate 
origin  is  the  blade  leading  edge.  The  following  quantities  are 
given:  blade  chord;  stagger  angle;  camber  angle;  section  area; 
location  of  centroid  of  the  section;  second  moments  of  area  of 
the  section  about  the  centroid;  orientation  of  the  principal 
axes;  and  the  principal  second  moments  of  area  of  the  section 
about  the  centroid.  Then  are  listed  the  coordinates  of  the 
camber  line,  the  camber  line  angle,  the  section  thickness,  and 
the  coordinates  of  the  blade  surfaces.  NPOINT  values  are  given. 
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A  lineprinter  plot  of  the  normalized  section  follows.  The  scales 
for  the  plot  are  arranged  so  that  the  section  just  fills  the  page 
so  that  the  scales  will  generally  differ  from  one  plot  to  another 
"Dimensional"  details  of  the  blade  section  are  given  next.  The 
normalized  data  given  previously  is  scaled  to  give  a  blade 
section  as  defined  by  IPCORD  and  CORD.  For  this  section  of  the 
output,  the  coordinates  are  with  respect  to  the  blade  stacking 
axis.  The  following  quantities  are  given:  blade  chord;  radius 
and  location  of  center  of  leading  (and  trailing)  edge(s); 
section  area;  the  second  moments  of  area  of  the  section  about  the 
centroid;  and  the  principal  second  moments  of  area  of  the 
section  about  the  centroid.  The  coordinates  of  NPOINT  points 
on  the  blade  surfaces  are  then  listed,  followed  by  the 
coordinates  of  31  points  distributed  at  (roughly)  six  degree 
intervals  around  the  leading  (and  trailing)  edges.  Finally, 
the  coordinates  of  the  blade  surfaces  and  points  around  the 
leading  (and  trailing)  edge(s)  is  (are)  shown  in  Cartesian 
form. 


The  quantities  required  for  aerodynamic  analysis  are 
printed  at  all  computing  stations  specified  by  the  IFANGS 
parameter.  The  radius,  section  angle,  blade  lean  angle,  blade 
blockage,  and  relative  angular  location  of  the  camber  line  are 
printed  at  each  streamsurface  intersection  with  the  particular 
computing  station. 

Details  of  the  manufacturing  sections  are  printed  if 
IPRINT  =  0  or  2.  At  each  value  of  z  specified  by  ZINNER, 

ZOUTER  and  NZ,  section  properties  and  coordinates  are  given. 

The  origin  for  the  coordinates  is  the  blade  stacking  axis. 

The  following  quantities  are  given:  section  area;  the  location 
of  the  centroid  of  the  section;  the  second  moments  of  area  of 
the  section  about  the  centroid;  the  principal  second  moments  of 
area  of  the  section  about  the  centroid;  the  orientation  of  the 
principal  axes;  and  the  section  torsional  constant.  Then 
the  coordinates  of  NPOINT  points  on  the  blade  section  surfaces 
are  listed,  followed  by  31  points  around  the  leading  (and 
trailing)  edge ( s ) . 

Precision  plots  are  produced  if  IFPLOT  =  1,  2,  3  or  *)  as 
described  under  the  definition  of  IFPLOT  given  previously. 

If  IPUNCH  =  1,  the  program  punches  the  quantities  required 
for  aerodynamic  analysis,  together  with  identifying  indices 
denoting  station  number  and  streamsurface  number,  on  cards  in 
the  following  format:  5  fields  each  of  12  locations  for  the 
quantities  themselves,  followed  by  2  fields  each  of  3  locations 
for  the  indices. 


SECTION  VIII 


EXAMPLES  OF  USE  OF  THE  PROGRAM 


1 .  INTRODUCTION 

This  section  shows  the  use  of  the  program  to  generate  four 
blades:  one  based  upon  the  polynomial  camber  line;  another 
upon  the  exponential  camber  line;  a  double-circular-arc  blade; 
and  one  based  upon  the  multiple-circular-arc  camber  line. 

The  same  basic  aerodynamic  design  is  assumed  for  the  first  two 
and  last  cases,  although  all  blades  could  not  satisfy  the  same 
aerodynamic  requirements.  Figure  18  is  a  meridional  section 
through  the  compressor  annulus,  showing  eight  streamsurfaces 
(including  the  hub  and  casing)  defining  the  flow,  and  the  locations 
of  six  planes  on  which  manufacturing  sections  are  to  be  determined. 
Usually,  more  streamlines  would  be  used  to  define  the  flow  and 
thus  the  blade,  by  means  of  the  streamsurface  blade  sections. 

Also,  it  would  normally  be  advantageous  to  determine  more  than 
six  manufacturing  sections.  However,  these  numbers  were  selected 
for  the  sample  runs  of  the  program  as  they  exemplify  the  program 
adequately  and  lend  clarity  to  the  resulting  figures.  Included 
on  Figure  18  are  the  approximate  locations  of  the  blade  leading 
and  trailing  edges  for  the  polynomial  camber  line  case.  The 
exponential  and  multiple-circular-arc  camber  line  cases  would 
appear  slightly  different  for  reasons  discussed  below. 

2.  POLYNOMIAL  CAMBER  LINE 

a.  Input  Data 

The  input  data  deck  used  for  this  example  is  listed 
below.  Some  points  of  interest  are  noted  in  the  order  in  which 
they  occur  in  the  input. 

As  mentioned  previously,  eight  streamsurface  blade 
sections  are  used  to  define  the  blade.  The  streamsurface  radii 
are  specified  at  eight  locations,  the  first  and  last  of  which 
are  outside,  and  well  clear  of,  the  rotor.  This  ensures  that 
the  boundary  conditions  imposed  on  the  spline-curve  (zero 
curvature  at  the  end  points)  have  little  influence  on  the  shape 
of  the  curve  representing  the  streamsurface  in  the  region  of 
interest,  that  is,  within  the  blade.  The  parameters  defining 
the  streamsurface  blade  sections  are  given  at  eight  locations; 
that  is,  directly  at  each  streamsurface.  Thus  the  interpolation 
of  the  parameters  at  each  streamsurface  is  reduced  in  this 
case  to  merely  reading  the  appropriate  row  of  data  from  the 
table.  Fifty  points  are  specified  to  define  each  blade  surface. 
This  number  serves  the  purposes  of  this  example  well,  but  it 
would  be  advantageous  to  use  more  points  where  the  precision- 
plot  output  is  to  be  directly  incorporated  in  the  manufacturing 
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procedure.  Specification  of  the  meridionally -projected  chord 
(rather  than  actual  chord)  is  selected.  All  optional  sections 
of  the  output  are  to  be  printed.  Superimposed  plots  of  the 
streainsurface  and  manufacturing  sections  are  to  be  produced  . at 
a  scale  of  five  times  full  size. 

Blade  inlet  and  outlet  angles  would  normally  be 
determined  from  an  aerodynamic  analysis  of  the  flow  through  trie 
blade  row.  Desirable  incidence  angles  for  this  profile  type 
are  probably  near  to  one  half  the  leading  edge  wedge  angle. 
Considerations  of  static  pressure  (or  velocity)  distributions 
within  the  blade  channel,  and  choking  flow  limits,  will  influenc 
the  selection  of  incidence  angles,  depending  upon  the  selected 
aerodynamic  analysis  method.  Deviation  angles  for  the  profile 
have  not  been  empirically  determined  as  yet,  and  therefore  the 
best  means  of  estimating  deviation  angles  is  probably  to  use 
Carter's  Rule,  or  a  similar  formulation.  Some  consideration  of 
the  camber  line  shape  (point  of  maximum  camber,  for  example) 
should  be  made.  The  parameters  PP  and  QQ  are  set  to  0.0  and  0.5 
on  all  s treamsur faces ,  respectively.  It  is  anticipated  that  the 
zero  curvature  leading  edge  configuration  (PP  =  0.0)  will  be 
frequently  specified  for  this  profile  type.  The  specification 
of  a  second  derivative  at  the  trailing  edge  that  is  less  than 
its  maximum  value  (QQ  =  0.5)  is  made  to  prevent  the  curvature 
of  the  camber  line  from  rising  to  a  maximum  at  the  trailing 
edge.  It  is  believed  that  this  could  cause  excessive  deviation 
angles.  The  blade  thickness  distribution  of  a  rotor  blade  is 
determined  by  aerodynamic  and  mechanical  factors.  The  leading 
edge  radius  and  trailing  edge  half-thickness  are  set  to 
approximately  .005  inch,  which  is  probably  a  practical  minimum. 
Maximum  thickness/chord  ratios  specified  vary  from  ^ . 8%  at  the 
hub  to  3.1$  at  the  casing.  These  figures  are  typical  of  those 
aerodynamically  acceptable  for  a  high  Mach  number  rotor  blade, 
and  yield  a  blade  with  a  reasonable  distribution  of  cross- 
sectional  area.  The  maximum  thickness  is  placed  at  0.6  of  the 
camber  line  from  the  leading  edge,  which  helps  to  maintain  a 
small  leading  edge  wedge  angle. 
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b .  Output  Data 

The  input  data  specified  all  optional  output  data, 
and  a  selection  for  each  section  of  the  output  is  shown  below. 
Some  points  of  interest  are  as  follows . 

Shown  first  is  the  printout  of  the  input  data;  all 
input  items  are  listed.  This  is  followed  by  details  of  the 
first  streamsurface  blade  section.  The  parameters  defining  the 
blade  section  are  followed  by  details  of,  and  a  line-printer 
plot  of,  the  normalized  blade  section.  Then  appears  a 
specification  of  the  section  scaled  to  the  desired  dimensions. 
First  there  are  the  fifty  points  specified  for  each  blade 
surface.  These  are  printed  two  per  line,  and  if  an  odd  number 
of  points  is  specified,  coordinates  of  the  last  point  are  not 
printed.  Next  follow  the  coordinates  of  the  31  points  describing 
the  leading  edge  radius.  The  final  data  for  the  streamsurface 
section  are  the  equivalent  Cartesian  coordinates  for  the  points 
just  mentioned.  It  will  be  seen  that  the  first  and  last  points 
on  the  leading  edge  radius  coincide  with  the  first  points  on 
the  pressure  and  suction  surfaces. 

The  same  format  is  repeated  for  each  of  the  remaining 
streamsurface  blade  sections,  and  a  printout  of  the  volume 
enclosed  by  the  blade,  but  these  results  have  not  been 
reproduced  here.  The  calculations  for  aerodynamic  analysis  are 
printed  next  for  those  computing  stations  for  which  they  were 
requested  by  the  1FANGS  parameter,  presented  in  order  of 
increasing  radial  streamsurface  location. 

Details  of  the  six  manufacturing  sections  defining 
the  blade  follow.  Reproduced  below  is  the  output  relating  to 
the  first  (innermost)  section.  Properties  of  the  section 
are  followed  by  the  coordinates  of  50  points  on  each  blade 
surface,  and  31  points  around  the  leading  edge.  It  will  be 
noted  that  the  section  centroid  is  not  calculated  to  be 
exactly  on  the  stacking  axis  (the  origin  for  the  coordinates) 
although  the  offsets  DELX  and  DELY  were  all  zero.  This  is  a 
reflection  of  the  fact  that  the  streamsurface  sections  is 
not  planar.  If  it  were  desired  to  have  the  centroids  of  the 
manufacturing  sections  lie  precisely  on  the  stacking  axis, 
the  program  would  be  rerun  with  DELX  and  DELY  offsets 
specified  so  that  they  counteract  the  mislocation  of  the 
centroids  previously  determined.  This  iterative  procedure 
might  require  more  than  one  loop,  depending  upon  the  accuracy 
desired.  Use  of  the  centroid  offsets  alo  permits  the  final 
blade  to  be  leaned  for  stress  redistribution  or  other  purposes. 

The  input  data  specified  superimposed  precision  plots 
of  both  the  streamsurface  sections  and  manufacturing  sections. 
These  plots  are  reproduced  (at  reduced  size)  as  Figures  19  and 
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20,  respectively.  It  is  of  interest  to  refer  also  to  Figure  18. 
The  innermost  manufacturing  plane  is  well  below  the  lowest  point 
on  the  hub  streamsurface  section  in  the  plane  of  the  stacking 
axis,  that  is,  the  plane  of  Figure  18.  The  Cartesian  coordinates 
of  the  first  streamsurface  blade  section  show  that  the  lowest 
point  on  the  section  (the  20th  point  on  the  leading  edge)  is  at 
z  =  6.66H 78.  The  radius  of  the  streamsurface  at  this  point  is 
6.75  (approximately),  and  the  innermost  manufacturing  plane  is 
at  z  =  6.5.  Thus,  at  the  leading  edge  the  extrapolation 
required  to  define  the  manufacturing  section  is  somewhat  smaller 
than  might  at  first  appear.  At  the  trailing  edge,  extrapolation 
is  required  to  define  the  first  two  manufacturing  sections. 

The  streamsurface  radius  at  the  casing  and  z-coordinate  for 
the  outermost  manufacturing  plane  both  equal  9.0  at  the  blade 
inlet.  However,  the  z-coordinate  of  the  leading  edge  of  the 
blade  is  8.8357,  so  that  the  outermost  manufacturing  section  too 
is  defined  completely  by  extrapolation.  Of  course,  portions 
of  the  blade  that  are  defined  by  extrapolation  would  not 
appear  on  a  final  blade,  but  would  facilitate  manufacture. 
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3.  EXPONENTIAL  CAMBER  LINE 


a.  Input  Data 

As  mentioned  previously,  the  overall  design  incorporated 
in  the  polynomial  camber  line  example  was  also  assumed  for  the 
exponential  camber  line  example.  The  discussion  of  the  input 
data  therefore  highlights  the  changes  made  to  the  polynomial 
camber  line  example  data  in  order  to  create  the  exponential  camber 
line  example. 

The  exponential  camber  line  is  specified  by  setting 
ISECN  equal  to  1.  For  this  example,  chords  rather  than 
meridionally-projected  chords  are  specified.  The  chord  lengths 
used  are  taken  from  the  results  of  the  polynomial  camber  line 
example,  and  because  of  the  different  camber  line  configurations, 
the  resulting  meridionally-projected  chords  are  somewhat 
smaller  for  the  exponential  camber  line  case.  The  parameters  PP 
and  QQ  are  set  equal  to  0.5  for  all  sections.  These  values 
will  probably  be  satisfactory  for  most  applications  of  the 
blade  profile.  The  only  remaining  change  to  the  polynomial 
camber  line  input  data  is  to  specify  the  distributions  of  the 
inflection  point  and  change  in  camber  angle  from  blade  inlet 
to  the  inflection  point.  The  distributions  used  for  the 
example  are  quite  arbitrary,  but  representative  of  what  might 
be  indicated  by  a  typical  aerodynamic  analysis.  The  inflection 
point  is  located  at  the  leading  edge  on  the  innermost  stream- 
surface  section,  and  moved  rearward  .05  per  streamsurface ,  to 
•35  (of  the  meridionally-projected  chord)  at  the  casing. 

(As  a  value  of  s  of  zero  is  not  acceptable,  the  value  used  for 
the  innermost  section  is  .01.)  For  the  three  innermost 
sections,  the  inflection  angle  is  set  equal  to  the  inlet  angle. 
Thus,  these  sections  have  straight  leading  edge  portions  (of 
length  from  zero  to  .1  of  the  meridionally-projected  chord), 
followed  by  portions  of  continuously-positive  camber.  The 
remaining  sections  have  inflection  angles  that  are  greater  than 
the  section  inlet  angles,  the  difference  in  angles  rising  from 
zero  at  the  third  section  (as  described  above)  to  seven  degrees 
at  the  casing.  Thus,  these  sections  have  reverse  and  positive 
camber,  the  reverse  camber  increasing  with  radius,  as  the 
inflection  point  also  moves  rearward. 

The  input  data  deck  is  listed  below. 
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b.  Output  Data 

The  input  calls  for  all  optional  printed  output.  As 
this  is  identical  in  format  to  that  shown  for  the  polynomial 
camber  line  example,  none  has  been  reproduced  here. 

Figures  21  and  22  show  (at  reduced  size)  the  super¬ 
imposed  streamsurface  section  and  manufacturing  section  plots. 

The  comments  regarding  extrapolation  made  in  connection  with 
the  polynomial  camber  line  example  again  apply.  A  comparison 
of  some  intermediate  sections  from  both  plots  shows  well  the 
value  of  an  accurate  method  of  determining  manufacturing 
coordinate  data  from  streamsurface  design  data.  For  blades 
where  the  streamsurface  section  form  is  critical  throughout 
the  section,  any  approximate  method  based  upon  blade  edge 
geometry  only  is  unlikely  to  be  satisfactory. 

4 .  DOUBLE-CIRCULAR-ARC  BLADE 

a.  Input  Data 

This  example  illustrates  several  program  features  not 
illustrated  by  the  previous  examples.  First,  this  example 
illustrates  curvilinear  computing  stations.  Secondly,  the 
blade  is  described  in  terms  of  negative  inlet  angle,  illustrating 
that  the  program  accommodates  blades  specified  in  either  sense. 
This  allows  stators  to  be  specified  in  the  opposite  sense  to 
rotors,  if  desired.  Finally,  the  blade  is  stacked  at  its 
trailing  edge,  rather  than  at,  or  offset  slightly  from,  the 
centroid. 


The  double-circular-arc  blade  is  specified  by  setting 
ISECN  equal  to  2.  The  particular  example  illustrates  a  novel 
design  concept  in  that  the  leading  edge  is  swept  from  hub  and 
tip  toward  mid-passage,  a  concept  which  is  described  in  further 
detail  in  Reference  4. 

The  input  data  is  listed  below. 
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b .  Output  Data 

The  printed  output  for  this  example  is  similar  in  most 
respects  to  that  shown  for  the  polynomial  camber  line.  The 
section  entitled  "Dimensional  Results"  is  shown  for  stream- 
surface  section  1  to  illustrate  that  the  section  coordinates 
are  presented  relative  to  the  stack  axis  which,  in  this  case, 
has  been  specified  at  the  blade  trailing  edge.  Two  other  minor 
differences  in  comparison  to  the  polynomial  camber  line  output 
are  related  to  the  trailing  edge:  The  printout  contains  the 
coordinate  location  of  the  center  of  the  trailing  edge  as  well 
as  the  31  points  distributed  around  the  trailing  edge.  The 
latter  feature  also  appears  with  the  printout  of  the  manu¬ 
facturing  sections  (not  reproduced  here). 

Figures  23  and  2k  show  (at  reduced  size)  streamsurface 
section  and  manufacturing  section  plots,  respectively.  The 
comments  regarding  extrapolation  made  in  connection  with  the 
polynomial  camber  line  example  again  apply. 
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POINTS  DESCRIBING  LEADING  ANC  TRAILING  ECGES 
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CARTESIAN  CCCRCINATES  ON  STREAMSLRFACE 
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5. 


MULTIPLE-CIRCULAR-ARC  CAMBER  LINE 


a.  Input  Data 

The  same  overall  design  incorporated  in  the  first 
two  examples  was  also  assumed  for  the  multiple-circular-arc 
camber  line  example.  The  discussion  of  the  input  data  will  be 
limited  to  changes  made  to  the  exponential  camber  line  example 
to  create  the  multiple-circular-arc  camber  line  example. 

The  multiple-circular-arc  camber  line  is  specified  by 
setting  ISECN  equal  to  3.  The  only  other  changes  made  to  the 
exponential  camber  line  example  are  the  distributions  of 
inflection  point  and  the  change  in  camber  angle  from  blade 
inlet  to  the  inflection  point.  The  inflection  point  has  been 
arbitrarily  set  at  0.3  (of  the  meridionally-projected  chord) 
on  all  streamsurfaces .  The  change  in  camber  angle  has  been 
specified  in  such  a  way  that  the  blade  has  positive  camber  in 
the  leading  segment  at  the  hub,  gradually  decreasing  to  zero 
camber  (straight  line  leading  segment)  on  streamsurface  *),  and 
producing  sections  with  reverse  camber  on  the  four  outermost 
streamsurfaces . 

The  input  data  deck  is  listed  below. 


FIGURE  MULTIPLE  CIRCULAR  ARC  CAMBER  LINE 
8868  50  30  203130 
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6.75 
7.0685 
7  .3844 
7.6994 
8.0155 
8.3365 
8.6625 
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1  1 

-.5767  0.0 
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e.0063 
8.2777 
8  »555f 
8.8417 
1  1 
.20003 
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7.5164 
7.7436 
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7.5766 


0.G 
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7 .736  C 
7.9054 
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b .  Output  Data 

The  input  data  calls  for  all  optional  printed  output. 
As  this  is  identical  in  format  to  that  shown  for  the  polynomial 
camber  line,  none  has  been  reproduced  here. 

Figures  25  and  2 6  show  (at  reduced  size)  the  super¬ 
imposed  streamsurface  section  and  manufacturing  section  plots, 
respectively.  The  comments  regarding  extrapolation  made  in 
connection  with  the  polynomial  camber  line  example  are  again 
applicable . 
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SECTION  IX 


COMPUTER  PROGRAM  DETAILS 


1.  IMPLEMENTATION  OF  THE  COMPUTER  PROGRAM 

The  program  is  written  in  FORTRAN  IV  and  was  developed  on 
an  IBM  709 V70^  Direct  Couple  System  (incorporating  the 
version  13  IBSYS  Operating  System)  and  a  CDC  6000  Series 
System  (incorporating  the  version  3-3  SCOPE  Operating  System). 
When  loaded  into  core,  the  program  (and  resident  system) 
occupies  about  32K  of  storage,  so  that  the  program  will  probably 
not  be  usable  without  modification  on  a  relatively  small 
computer.  Apart  from  this  limitation,  the  program  should  be 
compatible  with  the  majority  of  modern  computing  systems.  The 
program  consists  of  a  main  program,  and  Subroutines  BQ,  CQ,  Dl, 
EQ,  FQ  and  GQ.  The  seven  decks  have  been  given  the  identifiers 
A,  B,  C,  D,  E,  F  and  G,  respectively.  Listings  of  the  decks 
are  shown  below,  and  the  deck  set-up  required  for  the  CDC  6000 
Series  System  is  also  presented. 

The  program  uses  three  numerical  system  units  for  its 
input  and  output  routines.  Input  is  drawn  from  the  card  unit 
via  READ  statements  referring  to  Unit  LOG  1;  output  is  sent  to 
the  line  printer  by  WRITE  statements  referring  to  Unit  LOG  2; 
and  punched  output  is  produced  via  WRITE  statements  referenced 
to  Unit  LOG  3.  Units  LOG  1,  LOG  2,  and  LOG  3  are  set  equal  to 
5,  6,  and  7,  respectively  on  cards  All^O-60  in  the  FORTRAN 
programming.  On  the  CDC  6000  Series  System,  the  "PROGRAM"  card 
must  also  establish  the  input  and  output  linkages.  On  other 
computing  systems,  the  "PROGRAM"  card  may  not  be  required,  and 
the  input-output  files  may  be  established  via  control  cards. 

The  program  as  presented  herein  utilizes  an  on-line 
precision  plotting  capability  available  at  the  program 
development  site.  Calls  to  four  subroutines  not  included  in 
the  deck  are  included  in  the  program.  These  calls  are  executed 
only  if  precision  plots  are  specified  in  the  input  data,  and  the 
entry  points  expressed  are  part  of  standard  CALCOMP  software 
normally  supplied  to  users  of  CALCOMP  precision  plotting 
equipment.  The  various  call  statements  used  in  the  program 
are  explained  below,  to  facilitate  modification  should  the  need 
arise . 

CALL  PLOT  ( XPLOT,  YPLOT,  N) 

The  majority  of  plotting  is  done  using  this  form  of  call.  The 
parameters  XPLOT  and  YPLOT  are  the  "x"  and  "y"  coordinates  (in 
inches)  on  the  paper  to  which  the  pen  is  being  directed.  The 
parameter  N  indicates  pen  up  or  down,  |N|  =  3  or  |  N  j  =  2 
respectively,  and  will  cause  XPLOT  or  YPLOT  to  be  assigned 
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as  the  origin  for  further  coordinates  if  N  is  negative. 

CALL  SYMBOL  (X,  Y ,  H,  TEXT,  THETA,  N) 

This  call  is  used  to  title  the  plots.  The  parameters  X  and  Y 
are  the  coordinates  (in  inches)  of  the  lower  left  hand  corner 
of  the  first  character,  H  is  the  character  height  (in  inches), 
TEXT  is  the  character  to  be  printed,  THETA  is  the  angle  of  the 
lettering  with  respect  to  the  "x"  axis  and  N  is  the  total 
number  of  characters  to  be  printed. 

CALL  NUMBER  (X,  Y,  H,  P,  THETA,  N) 

This  call  causes  the  printing  of  the  number  P.  The  parameters 
X,  Y,  H,  and  THETA  are  used  as  for  CALL  SYMBOL.  The  parameter 
N  indicates  the  number  of  digits  following  the  decimal  point 
if  positive,  or  truncation  to  an  integer  if  equal  to  -1. 

CALL  PLOTE 

This  call  terminates  the  tape. 

In  the  event  the  program  is  used  on  a  computing  system 
which  does  not  include  CALCOMP  software,  and  the  operating 
system  will  not  execute  a  program  with  unsatisfied  external 
references,  dummy  entry  points  may  be  supplied  by  adding  to 
the  deck  a  subroutine  such  as  the  following: 

SUBROUTINE  PLOT 

A  =  A 

ENTRY  SYMBOL 
ENTRY  NUMBER 
ENTRY  PLOTE 
RETURN 
END 

2.  DECK  SETUP  FOR  CDC  6000  SERIES  SYSTEM 

The  deck  setup  required  to  run  the  program  on  a  CDC  6000 
Series  System  incorporating  the  SCOPE  3.3  Operating  System  is 
shewn  below.  Production  runs  of  the  program  would  usually 
employ  relocatable  binary  forms  of  the  routines  produced 
from  the  source  decks  to  avoid  having  to  compile  the  FORTRAN 
for  each  run  and  the  waste  of  associated  computer  resources. 
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JOB 


Identification,  etc. 


FTN. 

LGO. 

7/8/9  End  of  Record 
SOURCE  DECK  A 
SOURCE  DECK  B 
SOURCE  DECK  C 
SOURCE  DECK  D 
SOURCE  DECK  E 
SOURCE  DECK  F 


SOURCE  DECK  G 
7/8/9  End  of  Record 
DATA  DECK 

6/ 7/8/9  End  of  Job 
3.  FORTRAN  PROGRAM  LISTING 
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560  CALL  01  (SS,X,1uj,SS1(1,1> ,SSi(l,i)  ,N  ,1)  B  6960 

CALL  Dl (SS,X, 10u,SSi(l, 2) ,SSl(l,2) ,N  ,1)  B  6970 

CALL  01 (SS,X,lju,SSl(i,3) ,SS1(1,3> ,31  ,1)  B  6990 

IF(ISTAK.GT.l)  GO  TO  57C  B  6990 

IF(ISTAK.EQ.i)  SIGMAO=SIGMA (100)  B  7000 
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The  analysis  which  has  been  described  in  this  report  is 
performed  primarily  in  the  main  program  and  Subroutine  BQ, 
while  the  other  subroutine  performs  specific  tasks  supportive 
of  the  desired  objectives.  Subroutine  BQ  concerns  the  details 
of  defining  and  describing  the  blade  on  the  streamsurface;  the 
main  program  utilizes  this  information  to  stack  the  blade  and 
obtain  the  manufacturing  section  description.  Subroutine  CQ 
is  used  to  find  slopes  of  various  ’’spline-curves"  at  particular 
points.  D1  is  the  curve-fitting  routine,  EQ  produces  a  plot 
in  the  printed  portion  of  the  output,  and  PQ  is  used  to  draw 
streamsurface  and/or  manufacturing  section  plot  axes  if  IFPLOT  = 

1,  2  or  3.  Subroutine  GQ  is  used  to  compute  the  center  coordinate, 
the  radius  of  curvature,  and  the  arc  length  for  the  circular- 
arc  and  multiple-circular-arc  camber  lines. 

A  description  of  the  calculation  procedure  employed  in 
the  main  program  and  in  Subroutine  BQ  is  described  below. 

The  steps  of  this  procedure  are  keyed  to  their  location  in  the 
program  through  the  associated  deck  serialization,  which  appears 
parenthetically . 

1.  The  input  data  is  read  and  printed.  (A1190-A2160) 

2.  If  precision  plotting  is  specified,  the  plot  is  initialized 
and  if  IFPLOT  =  1  or  3,  the  axes  for  the  superimposed  stream- 
surface  plots  are  produced.  (A21T0-A2200 ) 

3.  A  loop  performed  for  each  streamsurface  section  is  commenced. 
This  loop  creates  the  blade  on  the  streamsurface.  (A2210) 

If.  The  axial  locations  of  the  intersections  of  a  particular 
streamsurface  with  all  computing  stations  are  determined,  the 
interval  between  the  first  and  last  points  is  subdivided  into 
99  intervals,  and  the  axial  coordinates  of  the  100  resulting 
locations  obtained.  ( A2220-A2290 ) 

5.  The  slope  of  the  streamsurface  at  each  of  the  100  locations 
obtained  in  If.  above  is  derived.  (A2300) 

6.  The  slope  of  the  computing  stations  at  each  of  the  stream- 
surface  computing  station  intersection  points  is  derived. 

(A2310) 

7.  Using  Equation  (7*0,  the  streamsurface  length  (in  the 
meridional  projection)  is  obtained  over  the  99  intervals  of 

If.,  providing  an  x-m  table  for  the  streamsurface.  (A2320-A23*J0 ) 

8.  The  parameters  defining  the  streamsurface  blade  section  are 
Interpolated  (or  extrapolated)  from  the  input  tables.  If 
NSPEC  =  1,  they  are  taken  to  be  radially  uniform.  If  NSPEC  =  2, 
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linear  interpolation  is  employed.  If  NSPEC  >2,  spline-curve 
interpolation  is  used.  (A2350-A2500) 

9.  The  design  of  the  streamsurface  section  is  initiated. 
(A25'10-A25*10) 

10.  If  IPRINT  =  0  or  1,  the  parameters  defining  the  stream- 
surface  blade  section  are  printed.  (131230-131*110) 

11.  If  ISECN  =  0,  the  coefficients  of  the  quartic  camber  line 
are  computed  using  Equations  (6),  (7),  (8),  (9)  and  (10). 

The  chord  length  corresponding  to  a  meridional  camber  line 
chord  of  unity  is  determined.  (B1*J50-B1510) 

12.  If  ISECN  =  1,  the  coefficients  for  the  two  segments  of  the 
exponential  camber  line  are  computed  using  Equations  (17),  (18), 
(19),  (20),  (21),  (22),  (23)  and  (2*1).  Equations  (18)  and  (2*1) 
require  iterative  solution,  and  up  to  50  attempts  to  satisfy 
the  equations  are  permitted.  The  equations  are  taken  to  be 
satisfied  when  the  parameter  P  or  Q,  as  appropriate,  is  given 

to  within  0.0001  of  its  specified  value.  (Should  the  iterative 
procedure  fail,  a  diagnostic  is  printed  and  the  job  is  terminated. 
This  will  not  normally  occur.)  The  chord  length  corresponding 
to  a  meridional  camber  line  chord  of  unity  is  determined. 
(B1530-B2050) 

13.  If  ISECN  =  3,  the  constants  of  the  equations  for  the  two 
segments  of  the  multiple-circular-arc  camber  line  are  determined 
from  Equations  (38),  (36),  (35),  and  (3*1)  respectively,  in 
Subroutine  GQ.  These  constants  are  scaled  and  the  arc  center 
offset,  if  required,  and  the  arc  length  of  each  segment  determined 
cumulatively.  The  chord  length  corresponding  to  a  meridional 
camber  line  chord  of  unity  is  determined.  (B2070-B2300) 

NOTE:  Each  segment  is  described  by  a  proportion  of  the  NPOINT 
camber  line  points  based  on  its  proportion  of  the  meridionally 
projected  chord  length.  If  this  number  of  points  for  the 
leading  segment  is  0  or  1,  this  segment  is  abandoned,  and 
one  circular  arc  is  produced  for  the  entire  camber  line. 

1*1.  If  ISECN  =  2,  the  constants  of  the  equation  for  the 
circular-arc  camber  line  and  the  camber  line  arc  length  are 
computed  as  for  ISECN  =  3.  The  chord  length  corresponding  to  a 
meridional  chord  of  unity,  that  corresponding  to  an  overall 
blade  chord  of  unity,  and  an  associated  scale  factor 
(described  in  15.)  are  determined  for  the  double-circular-arc 
blade.  (B2320-B2350) 

15.  For  ISECN  =  0,  1  or  3,  the  chord  length  corresponding  to 
an  overall  blade  chord  of  unity  as  well  as  a  scaling  factor 
which  relates  the  ratio  of  dimensions  for  an  overall  meridional 
chord  of  unity  to  dimensions  for  a  camber  line  meridional 


chord  of  unity  are  determined.  (B2370-B2380 ) 


16.  If  IFCORD  =  1,  the  specified  chord  is  divided  by  the  chord 
corresponding  to  a  unity  meridional  chord,  giving  the  desired 
meridional  chord.  The  three  blade  thickness  descriptors  are 
scaled  so  that  they  apply  to  a  blade  having  a  camber  line 
meridional  chord  of  unity.  (B2390-B2420) 

17*  The  coefficients  of  the  two  thickness  equations  of  the 
standard  thickness  distribution  are  computed  using  Equations 
(42),  (43),  (44),  (47),  (48)  and  (49).  (B2470-B2520) 

18.  The  length  of  the  section  camber  line  for  ISECN  =  0  or  1 
is  determined  by  numerical  integration.  The  meridional  chord 
(of  unity)  is  divided  into  (NPOINT-1)  uniform  intervals,  where 
NPOINT  is  the  number  of  points  specified  in  the  input  to  define 
each  blade  surface.  The  integration  for  the  curve  length 
between  each  adjacent  pair  of  points  is  obtained  using  Simpson’s 
Rule,  the  interval  being  subdivided  into  10  equal  intervals. 

If  the  polynomial  camber  line  is  specified,  the  gradient  of 
the  camber  line  is  given  by  Equation  (4).  If  the  exponential 
camber  line  is  specified,  the  gradient  of  the  camber  line  is 
given  by  Equation  (15),  different  coefficients  applying  forward 
and  rearward  of  the  inflection  point.  These  data  enable  a 
location  on  the  camber  line  to  be  expressed  as  a  fraction  of 
camber  line  length.  (B2540-B2730) 

19.  At  each  of  the  NPOINT  points,  equispaced  on  the  meridional 
chord,  the  section  half-thickness  is  computed  from  Equations 
(41),  (46)  or  (67),  and  the  coordinate  and  slope  of  the  camber 
line  are  computed,  from  either  Equations  (5)  and  (4)  (polynomial), 
(16)  and  (15)  (exponential),  or  (26)  and  the  derivative  of  (26) 
solved  for  (x-xQ),  choosing  the  appropriate  branch  of  the 

square  root  involved,  for  the  circular-arc  and  multiple-circular- 
arc  camber  lines.  (B2970-B3470 ) 

20.  The  coordinates  of  NPOINT  points  on  each  blade  surface  are 
determined  by  combining  the  data  derived  in  19.  The  scaling 
factor  derived  in  14 .  or  15.  is  applied  so  that  the  blade 
coordinates  correspond  to  an  overall  blade  meridional  chord  of 
unity.  (B3480-B3520) 

21.  The  slope  of  the  camber  line  is  expressed  in  degrees.  The 

camber  line  coordinates,  lengths,  and  the  section  half-thickness 
are  scaled  by  the  factor  derived  in  14.  or  15.  (B3530-B3570) 

22.  The  section  area  and  the  location  of  the  centroid  are 

determined  using  Equations  (8l),  (82)  and  (83).  (B3580-B3730) 

23.  The  camber  line  is  redefined  in  terms  of  100  points  to 
assure  sufficient  points  for  accurate  linear  interpolations 
in  the  determination  of  e,  needed  for  the  eventual  trans- 
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formation  onto  manufacturing  sections.  (63780-642*10) 

24.  The  scaled  camber  line  is  interpolated  to  yield  its 
description  over  99  equal  axial  intervals.  (B4250-B4260) 

25.  The  camber  line  y-coordinate  and  slope  at  the  axial 
location  of  the  centroid  or  centroid-offset  are  determined  by 
linear  interpolation.  (B4270-B4280) 

26.  The  product  of  inertia  and  second  moments  of  area  of  the 

section  about  the  centroid  are  determined  using  Equations  (8*0, 
(85),  (86),  (87),  (88),  (89),  (90)  and  (91).  (B4310-B4470) 

27.  The  orientation  of  the  principal  axes  and  the  principal 

second  moments  of  area  of  the  section  about  the  centroid  are 
determined  using  Equations  (92),  (93)  and  (9*0.  (B4480-R4520) 

28.  If  IPRINT  =  0  or  1,  details  of  the  normalized  blade  section 
(meridional  chord  unity)  are  printed.  (B4570-B4850) 

29.  If  IPRINT  =  0  or  1,  a  line-printer  plot  of  the  section  is 
made.  (B4870-B4940) 

30.  Section  properties  determined  in  22.,  26.,  and  27.  are 
scaled  according  to  the  specified  section  chord  to  produce 
"dimensional"  results.  (B4950-B5120) 

31.  If  IPRINT  =  0  or  1,  section  properties  are  printed. 
(B5150-B5*120) 

32.  The  coordinates  of  the  NP0INT  points  on  the  blade  surfaces 
are  scaled  according  to  the  specified  chord,  and  the  origin  is 
shifted  to  the  stacking  axis.  The  stacking  axis  is  offset 
from  the  centroid  of  the  section  by  the  specified  distances. 

If  IPRINT  =  0  or  1,  the  coordinates  are  orinted.  (B5510- 
B5670) 

33*  The  coordinates  of  31  points  describing  the  blade  edge(s) 
are  determined.  For  the  standard  thickness  distribution,  it  is 
assumed  that  the  wedge  angle  is  zero,  so  that  the  radius 
extends  over  a  complete  semicircle,  hence  the  points  are 
distributed  at  6°  intervals.  For  the  double-circular-arc 
blade,  the  2nd  and  30th  points  of  the  edges  are  determined  at 
the  exact  location  of  the  tangency  between  the  particular  blade 
surface  and  the  edge  radius,  resulting  in  an  angular  interval 
distribution  uniform  between  points  2  and  30,  with  irregular 
first  and  last  intervals.  (B5740-B6220) 

34.  The  origin  of  the  m-coordinate  system  is  shifted  to  the 
section  centroid  or  centroid-offset  as  provided  by  the  input 
data.  (B6250-B6290) 
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35.  The  m-length  corresponding  to  each  of  the  100  points 
describing  the  camber  line  is  determined.  (B6360-B6410) 

36.  The  reference  of  Equation  (77)  is  determined.  (B6420) 

37.  The  integration  of  Equation  (78)  is  performed  toward  the 
ends  of  the  blade  section  from  the  centroid  or  centroid-offset, 
giving  the  e  coordinate  for  the  camber  line.  (B6430-B6840) 

38.  The  m-coordinate  of  each  point  of  the  blade  surface  is 
determined.  (B687O-B698O ) 

39.  The  e  coordinate  system  is  shifted  to  the  leading  or 
trailing  edge  if  the  blade  is  stacked  there.  (B7000-B7030) 

*10.  The  e  corresponding  to  each  point  on  the  section  surface 
is  determined  from  the  e  of  the  camber  line  at  that  axial 
location  and  Ae  calculated  from  the  displacement  from  the 
camber  line  and  the  streamsurface  radius.  (B7070-B7260) 

41.  The  location  of  the  stacking  axis  is  given  (in  the  input 

data)  by  specifying  its  axial  coordinate.  Linear  interpolation 
yields  the  corresponding  streamsurface  coordinate  from  the  x-m 
table  created  in  7.  (A2550) 

42.  The  origin  of  the  m-coordinate  shifted  in  34.  is  again 

shifted,  to  the  stacking  axis  (determined  in  4l.).  (A2560- 

A2580) 

43.  The  origin  of  the  table  of  axial  coordinates  of  the 
streamsurface-computing  station  intersection  points  (obtained 
in  4.)  is  shifted  to  the  stacking  axis.  (A2590-A2600) 

44.  If  IPPLOT  =  1  or  3>  the  streamsurface  blade  section  plot 
is  produced.  (A2620-A2850) 

45.  A  trigger  is  set  if  the  calculation  of  quantities  for  aero¬ 
dynamic  analysis  of  the  blade  is  specified  at  any  computing 
station.  ( A2860-A2890 ) 

46.  If  the  calculations  for  aerodynamic  analysis  are  specified, 
the  angular  location  of  the  camber  line  with  respect  to  the 
stack  axis,  the  slope  of  the  camber  line  on  the  streamsurface, 
and  the  Cartesian  coordinates  of  the  camber  line  at  the  stream- 
surface-computing  station  intersections  are  determined  and 
stored.  ( A2910-A2960) 

47.  The  axial  coordinate  of  each  of  the  NPOINT  points  specifying 
the  streamsurface  blade  section  suction  surface  is  obtained 

by  linear  interpolation  from  the  x-m  table  produced  in  42. 
(A2970-A2990) 
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48.  The  streamsurface  radius  at  each  axial  location  derived 
in  47.  is  obtained  by  spline-curve  interpolation  from  the  table 
created  in  43*  (A3000) 


49.  Cartesian  coordinates  for  each  point  on  the  streamsurface 
section  suction  surface  are  obtained  using  the  ’  e's  computed 
in  40.  and  Equation  (79).  (A3010-A3060 ) 


50.  Steps  47.,  48.,  and  49.  are  repeated  for  the  pressure 
surface.  (A3070-A3160) 


51.  Steps  47.,  48.,  and  49.  are  repeated  for  the  points 
describing  the  leading  edge.  (A3170-A3260) 


52.  Steps  47.,  48.,  and  49.  are  repeated  for  the  points 
describing  the  trailing  edge  if  ISECN  =  2.  (A3280-A3370) 

53*  If  IPRINT  =  0  or  1,  the  Cartesian  coordinates  of  each 
point  describing  the  streamsurface  blade  section  are  printed. 
(A3390-A3910) 


54.  The  loop  initiated  in  3.  for  each  streamsurface  is 
terminated.  Upon  completion  of  this  loop,  the  Cartesian 
coordinates  of  all  streamsurface  sections  have  now  been  stored 
for  interpolation  on  the  specified  manufacturing  planes. 

(A 39 20) 


55.  Unless  IPRINT  =  1,  the  volume  of  the  blade  is  calculated, 
using  Equation  (106),  and  printed.  (A3950-A413O ) 

56.  If  the  calculations  for  aerodynamic  analysis  are  specified, 

the  program  prints  the  appropriate  heading.  If  not,  the 
program  moves  to  step  62.  (A4l60-A4l80) 

57.  In  the  original  x-r  coordinate  system,  the  axial  coordinate 
of  the  streamsurfaces-computing  station  intersections  and  the 
computing  station  slope,  are  determined  for  each  station  at 
which  the  aerodynamic  analysis  quantities  were  specified. 
(A4260-A4270) 

58.  The  blade  blockage  of  Equation  ( 10 7 )  is  obtained  at  each 
of  these  points.  (A4280-A4400) 

59.  The  blade  lean  angle  of  Equation  (108)  is  obtained  at  each 
of  these  points.  (A4510) 

60.  The  mean  section  angle  of  Equation  (109)  is  obtained  at 
each  of  these  points.  (A4530-A4540) 

61.  The  quantities  required  for  aerodynamic  analysis  are 

printed,  and  produced  on  cards  if  IPUNCH  =  1.  (A4550-A4570) 
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62.  If  IFPLOT  =  2  or  3,  the  axes  are  drawn  and  titled  for  the 
superimposed  plot  of  all  the  manufacturing  sections.  (A^620- 
A4630) 

63.  If  no  output  relating  to  manufacturing  sections  is 
specified  by  either  IFPLOT  or  IPRINT,  the  remainder  of  the 
program  is  bypassed.  Alternatively,  if  printed  details  of  the 
manufacturing  sections  are  specified,  a  heading  is  printed. 
(A*l61i0-A>4700) 

6k.  The  location  of  each  of  the  manufacturing  planes  is 
determined.  ( Ak710-Ak760 ) 

65.  The  (Cartesian)  coordinates  of  each  of  NPOINT  points  on 
the  blade  surface  are  obtained  by  spline-curve  interpolation 
at  each  of  the  manufacturing  sections.  ( A4770-Ai<860) 

66.  The  (Cartesian)  coordinates  of  each  of  31  points  describing 
the  blade  edge(s)  are  obtained  by  spline-curve  interpolation 

at  each  of  the  manufacturing  sections.  (Ail870-AiJ990) 

67.  A  loop  that  is  performed  for  each  manufacturing  section 
is  initiated.  This  loop  contains  the  determination  of  section 
properties  and  the  output  of  results  for  the  section.  (A5000) 

68.  The  section  area  and  location  of  the  centroid  are 

determined  using  Equations  (98),  (99)  and  (100).  (A5010-A5250) 

69.  The  torsional  constant  is  determined  using  Equations  (10*0 

and  (105).  (A5260-A5*I20) 

70.  The  second  moments  of  area  and  product  of  inertia  for  the 

manufacturing  section  about  its  centroid  are  determined  using 
Equations  (101),  (102),  (103),  (86),  (87),  (88),  (89),  (90) 
and  (91).  (A5^30-A5650) 

71.  From  the  second  moments  of  area  and  product  of  inertia, 

the  orientation  of  the  principal  axes  and  hence  principal 
moments  of  area  are  determined,  using  Equations  (92),  (93) 
and  (9*0  •  (A5660-A5690) 

72.  If  IPRINT  =  0  or  2,  section  properties  and  coordinates 
are  printed.  ( A5720-A6260 ) 

73-  If  IFPLOT  =  2  or  3,  a  plot  of  the  manufacturing  section  is 
produced.  ( A6310-A6530 ) 

7k.  If  IFPLOT  =  k ,  an  individual  plot  of  the  manufacturing 
sections  is  made.  The  axes  are  rotated  clockwise  until  the 
chord  line  is  horizontal.  The  angle  of  rotation  is  indicated 
as  the  stagger  angle.  (A6550-A7070) 
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75.  The  loop  initiated  in  Step  67.  for  each  manufacturing 
section  is  terminated.  (A7080) 

76.  If  precision  plots  have  been  made,  the  plotting  is 
terminated.  (A7090) 
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Figure  1  Variation  of  Normalized  Second  Derivative 

of  Polynomial  Camber  Line  with  Parameter  fP' 


L 


132 


\\\ 


•  • 


l-5-i 


Figure  5  Variation  of 
Alpha  2 

i: 


bu.  U 
0.00 
0.50 
10.0 
30.0 
50.0 


0.5  0.6  0.7  0.8  0.9  1.0 
X 


Polynomial  Camber  Line  with 


'/MOD (T ' * Mfl) 


60.00 

67.00 

40.00 

0.300 

0.300 

0.500 

0.700 


-0.2 


-0.4 


0.0  0.1  0.2  0.3  oTi  o75  oTi  0T7  578  0I9  T.o 

X 
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Figure  11  Some  Camber  Lines  Available  Under  the 

Multiple-Circular-Arc  Camber  Line  Option 


ALPHA  1 

ALPHAS 

ALPHA2 

S 

S 

S 


60.00 

67.00 

60.00 

0.300 

0.500 

0.700 


•  |f _  f  _ _ _ 

°Q  0.2  0.4  0.6  0.8  1.0 

X 


Figure  12  Variation  of  Multiple-Circular-Arc  Camber 
Line  with  Parameter  'S’ 
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Figure  13  Variation  of  Multiple-Circular-Arc  Camber 
Line  with  Parameter  'S' 
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Figure  14  Variation  of  Thickness  Distribution  with 
Point  of  Maximum  Thickness 
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STREAMSURFACE  BLADE  SECTION 


MANUFACTURING  BLADE  SECTION 


Figure  16  Calculation  of  Section  Properties 
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i\  AC 

(MBOL 

DESCRIPTION 

RELATED  INPUT  DATA  ITEM 

al 

Inlet  Angle 

B1 

a2 

Outlet  Angle 

B2 

*0 

Leading  Edge  Radius 

RLE  =  y0/c 

T 

Maximum  Thickness 

TC  =  T/C 

Trailing  Edge  Half -Thickness 

TE  =  yi/C 

Z 

Location  of  Maximum  Thickness 

Z  =  Z  (Camberline  Lenj 

C 

Chord 

CORD  if  IFCORD  =  1 

AC 

Axially-Projected  Chord 

CORD  if  IFCORD  =  0 

s 

Inflection  Point 

S  =  S/AC 

as 

Inflection  Point  Angle 

BS  +  B1 

Figure  IT  Parameters  Describing  Blade  Sections 
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Cartesian  Sections 


Figure  20  Example  Blade  Design  -  Polynomial 
Camber  Line 
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Streamsurface  Sections 


Figure  21  Example  Blade  Design 
Camber  Line 


Exponential 


Cartesian  Sections 


Figure  22  Example  Blade  Design  -  Exponential 
Camber  Line 
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Car*?  •■sian  Sections 
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Streamsurfaee  Sections 


Figure  25  Example  Blade  Design  -  Multiple- 
Circular-Arc  Camber  Line 
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Figure  26  Example  Blade  Design  -  Multiple- 
Circular-Arc  Camber  Line 
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